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This chapter discusses the equilibrium crystal shape (ECS) from a physical perspective, beginning 
with a historical introduction to the Wulff theorem. It takes advantage of excellent prior reviews, 
particularly in the late 1980s, recapping highlights from them. It contains many ideas and experi¬ 
ments subsequent to those reviews. Alternatives to Wulff constructions are presented. Controversies 
about the critical behavior near smooth edges on the ECS are recounted, including the eventual 
resolution. Particular attention is devoted to the origin of sharp edges on the ECS, to the impact 
of reconstructed or adsorbed surface phases coexisting with unadorned phases, and to the role and 
nature of possible attractive step-step interactions. 


Reformatted version (with slightly modified bibliog¬ 
raphy, also alphabetized and including article titles) 
of a review appearing in Handbook of Crys¬ 
tal Growth, Fundamentals, 2nd ed., edited 
by T. Nishinaga (Elsevier, Amsterdam, 2015-ISBN 
9780444563699/eBook:9780444593764), vol. lA (Ther¬ 
modynamics and Kinetics), chap. 5, pp. 215-264; 
http://www.sciencedirect.com/science/article/ 
pii/B9780444563699000058. 


1. INTRODUCTION 

The notion of equilibrium crystal shape is arguably 
the platonic ideal of crystal growth. It underpins much 
of our thinking about crystals and, accordingly, has 
been the subject of several special reviews and tutori¬ 
als [Rottman and Wortis 1984a, Wortis 1988, Zia 1988, 
Williams and Bartelt 1989, Bonzel 2003] and is a promi¬ 
nent section of most volumes and extended review articles 
and texts about crystals and their growth [Landau and 
Lifshitz 1980, van Beijeren and Nolden 1987, Nozieres 
1992, Pimpinelli and Villain 1998, Sekerka 2004]. In ac¬ 
tual situations, there are many complications that thwart 
observation of such behavior, including kinetic barriers, 
impurities, and other bulk defects like dislocations. Fur¬ 
thermore, the notion of a well-defined equilibrium shape 
requires that the crystal not make contact with a wall or 
surface, since that would alter its shape. By the same 
token, the crystal cannot then be supported, so gravity 
is neglected. For discussions of the effect of gravity or 
contact with walls, see, e.g., Nozieres [1992]. 

Gibbs [1874-8] is generally credited with being the 
first to recognize that the equilibrium shape of a sub¬ 
stance is that which, for a fixed volume, minimizes the 
(orientation-dependent) surface free energy integrated 
over the entire surface: the bulk free energy is irrelevant 
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since the volume is conserved, while edge or corner en¬ 
ergies are ignored as being higher-order effects that play 
no role in the thermodynamic limit. Herring [1951, 1953] 
surveys the early history in detail: The formulation of 
the problem was also carried out independently by Curie 
[1885]. The solution of this ECS problem, the celebrated 
Wulff construction, was stated by by Wulff (1901), but 
his proof was incorrect. Correct proofs were subsequently 
given by Hilton [1903], by Liebmann [1914], and by von 
Laue [1943], who presented a critical review. However, 
these proofs, while convincing of the theorem, were not 
general (and evidently applied only to T=0, since they 
assumed the ECS to be a polyhedron and compared the 
sum over the facets of the surface free energy of each facet 
times its area with a similar sum over a similar polyhe¬ 
dron with the same facet planes but slightly different ar¬ 
eas (and the same volume). Dinghas [1944] showed that 
the Brunn-Minkowski inequality could be used to prove 
directly that any shape differing from that resulting from 
the Wulff construction has a higher surface free energy. 
Although Dinghas again considered only a special class 
of polyhedral shapes. Herring [1951, 1953] completed the 
proof by noting that Dinghas’s method is easily extended 
to arbitrary shapes, since the inequality is true for con¬ 
vex bodies in general. In their seminal paper on crys¬ 
tal growth. Burton, Cabrera, and Frank [1951] present a 
novel proof of the theorem in two dimensions (2D). 

Since equilibrium implies minimum Helmholtz free en¬ 
ergy for a given volume and number, and since the bulk 
free energy is ipso facto independent of shape, the goal 
is to determine the shape that minimizes the integrated 
surface free energy of the crystal. The prescription takes 
the following form: One begins by creating a polar plot 
of the surface free energy as a function of orientation an¬ 
gle (of the surface normal) and draws a perpendicular 
plane (or line in 2D) through the tip of each ray. (There 
are many fine reviews of this subject by Chernov [1960], 
Frank [1962], Mullins [1962], and Jackson [1975].) Since 
the surface free energy in 3D (three dimensions) is fre¬ 
quently denoted by 7 , this is often called a 7 plot. The 
shape is then the formed by the interior envelope of these 
planes or lines, often referred to as a pedal. At zero tem¬ 
perature, when the free energy is just the energy, this 
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shape is a polyhedron in 3D and a polygon in 2D, each 
reflecting the symmetry of the underlying lattice. At fi¬ 
nite temperature the shapes become more complex. In 
2D the sharp corners are rounded. In 3D, the behavior 
is richer, with two possible modes of evolution with ris¬ 
ing temperature. For what Wortis terms type-A crystals, 
all sharp boundaries smoothen together, while in type- 
B, first the corners smooth, then above a temperature 
denoted Tq the edges also smooth. The smooth regions 
correspond to thermodynamic rough phases, with height- 
height correlation functions that diverge for large lateral 
separation I —like P, with a (typically 0 < a < 1 ) called 
the roughening exponent—in contrast to facets, where 
they attain some finite value as Z —>■ oo [Pimpinelli and 
Villain 1998]. The faceted regions in turn correspond 
to “frozen” regions. Pursuing the correspondence, sharp 
and smooth edges correspond to first-order and second- 
order phase transitions, respectively. 

The aim of this chapter is primarily to explore physical 
ideas regarding ECS and the underlying Wulff construc¬ 
tions. This topic has also attracted considerable interest 
in the mathematics community. Readers interested in 
more formal and sophisticated approaches are referred to 
two books, Dobrushin et al. [1992] and Cerf and Picard 
[2006] and to many articles, including De Coninck et al. 
[1989], Fonesca [1991], Pfister [1991], Fonesca and Muller 
[1992], Dacorogna and Pfister [1992], Dobrushin et al. 
[1993], Miracle-Sole and Ruiz [1994], Almgren and Tay¬ 
lor [1996], Peng et al. [1999], Miracle-Sole [1999, 2013]. 

2. FROM SURFACE FREE ENERGIES TO 
EQUILIBRIUM CRYSTAL SHAPE 

2.1. General Considerations 



FIG. 1: Portion of a (3,2,16) surface, vicinal to an fee (001), 
to illustrate a misoriented, vicinal surface. The vicinal-surface 
and terrace normals are n = (3, -2, 16)/\/269 and no = (0, 0, 
1), respectively. The polar angle 9 (with respect to the (001) 
direction), denoted cji in the original figure (consistent with 
most of the literature on vicinal surfaces), is arccos( 16/V269), 
while azimuthal angle (denoted 9 in most of the literature on 
vicinal surfaces), indicating how much h is rotated around ho 
away from the vertical border on which 9o is marked, is clearly 
arctan(l/5); tan^o = tan^ocosyi. Since h is where 

ai is the nearest-neighbor spacing, the mean distance I (in a 
terrace plane) between steps is ai/(%/2tan0) = 8-y/2/13ai = 
3.138ai. While the average distance from one step to the next 
along a principal, (110) direction looks like 3.5ai, it is in fact 
ai/(\/2 tan^o) = 3.2ai. The “projected area” of this surface 
segment, used to compute the surface free energy fp, is the 
size of a (001) layer: 20ai x 17ai = 340ai; the width is 20ai. 
In “Maryland notation” [see text] 2 is in the ho direction, 
while the normal to the vicinal, h, lies in the x — z plane and 
y runs along the mean direction of the edges of the steps. In 
most discussions, i/p = 0, so that this direction would be that 
of the upper and lower edges of the depicted surface. Adapted 
from Nelson et al. [1993] 


To examine this process more closely, we examine the 
free energy expansion for a vicinal surface, that is a sur¬ 
face misoriented by some angle 9 from a facet direction. 
Cf. Fig. 1. Unfortunately, this polar angle is denoted by cj) 
in much of the literature on vicinal surfaces, with 9 used 
for in-plane misorientation; most reviews of ECS use 9 
for the polar angle, as we shall here. The term “vicinal” 
implies that the surface is in the vicinity of the orien¬ 
tation. It is generally assumed that the surface orienta¬ 
tion itself is rough (while the facet direction is below its 
roughening temperature and so is smooth). We consider 
the projected surface free energy fp{9, T) [Jayaprakash et 
al. 1984] (with the projection being onto the low-index 
reference, facet direction of terraces) : 


/p (0, T) =/o (T) +/3 (T) + g(r) I tan 013 + c(tan 0)^ 

( 1 ) 

The first term is the surface free energy per area of the 
terrace (facet) orientation; it is often denoted a. The av¬ 
erage density of steps (the inverse of their mean separa¬ 
tion (£)) is tan0/h, where h is the step height. In the sec¬ 


ond term j3{T) is the line tension or free energy per length 
of step formation. (Since 2D is a dimension smaller than 
3D, one uses /3 rather than 7 . Skirting over the differ¬ 
ence in units resulting from the dimensional difference, 
many use 7 in both cases. While step free energy per 
length and line tension are equivalent for these systems, 
where the surface is at constant (zero) charge, they are 
inequivalent in electrochemical systems, where it is the 
electrode potential conjugate to the surface charge that is 
held fixed [Ibach and Schmickler, 2003]) The third term 
is associated with interactions between steps, in this case 
assumed to be proportional to (so that this term, 
which also includes the density of steps, goes like 
The final term is the leading correction. 

The interaction is due to a combination (not a sim¬ 
ple sum) of two repulsive potential energies: the entropic 
repulsion due to the forbidden crossing of steps and the 
elastic repulsions due to dipolar strains near each step. 
An explicit form for g{T) is given in Eq. (27) below. 
The of the entropic interaction can be understood 
from viewing the step as performing a random walk in 
the direction between steps (the x direction in “Mary- 
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land notation” [T] as a function of the y direction (which 
is time-like in the fermion transcription to be discussed 
later)—cf. Fig. 1, so the distance (y) it must go till it 
touches a neighboring step satishes oc y. To get a 
crude understanding of the origin of the elastic repulsion, 
one can imagine that since a step is unstable relative to 
a flat surface, it will try to relax toward a flatter shape, 
pushing atoms away from the location of the step by a 
distance decaying with distance from the step. When two 
steps are close to each other, such relaxation will be frus¬ 
trated since atoms on the terrace this pair of steps are 
pushed in opposite directions, so relax less than if the 
steps are widely separated, leading to a repulsive inter¬ 
action. Analyzed in more detail [Marchenko and Parshin 
1980, Nozieres 1992, Stewart et al. 1994], this repulsion 
is dipolar and so proportional to However, attempts 
to reconcile the prefactor with the elastic constants of 
the surface have met with limited success. The quartic 
term in Eq. Q is due to the leading correction to 

the elastic repulsion [Najafabadi and Srolovitz 1994], a 
dipole-quadrupole repulsion. It generally has no signif¬ 
icant consequences but is included to show the leading 
correction to the critical behavior near a smooth edge on 
the ECS, to be discussed below. 

The absence of a quadratic term in Eq. Q reflects that 
there is no interaction between steps. In fact, there 
are some rare geometries, notably vicinals to (110) sur¬ 
faces of fee crystals (Au in particular) that exhibit what 
amounts to repulsions that lead to more subtle be¬ 
havior [Carlon and van Beijeren 2000]. Details about this 
faxcinating idiosyncratic surface are beyond the scope of 
this chapter; readers should see the thorough, readable 
discussion by van Albada et al. [2002]. 

As temperature increases, /3(r) decreases due to in¬ 
creasing entropy associated with step-edge excitations 
(via the formation of kinks). Eventually /3(T) vanishes at 
a temperature Tr associated with the roughening tran¬ 
sition. At and above this Tr of the facet orientation, 
there is a profusion of steps, and the idea of a vicinal 
surface becomes meaningless. For rough surfaces the pro¬ 
jected surface free energy fp{9,T) is quadratic in tan0. 
To avoid the singularity at 0 = 0 in the free energy expan¬ 
sion that thwarts attempts to proceed analytically, some 
treatments, notably Bonzel and Preuss [1995], approxi¬ 
mate fp(0,T) as quadratic in a small region near 0 = 0. 
It is important to recognize that the vicinal orientation 
is thermodynamically rough, even though the underly¬ 
ing facet orientation is smooth. The two regions cor¬ 
respond to incommensurate and commensurate phases, 
respectively. Thus, in a rough region the mean spacing 
(£) between steps is not in general simply related to (i.e. 
an integer multiple plus some simple fraction) the atomic 
spacing. 

Details of the roughening process have been reviewed 
by Weeks [1980] and van Beijeren and Nolden [1987]; the 
chapter by Akutsu in this Handbook provides an up-to- 
date account. However, for use later, we note that much 
of our understanding of this process is rooted in the map¬ 


ping between the restricted BCSOS (body-centered [cu¬ 
bic] solid-on-solid) model and the exactly-solvable [Lieb 
1967, Lieb and Wu 1972] symmetric 6-vertex model [van 
Beijeren 1977], which has a transition in the same uni¬ 
versality class as roughening. (This BCSOS model is 
based on the BCC crystal structure, involving square net 
layers with ABAB stacking, so that sites in each layer 
are lateral displaced to lie over the centers of squares in 
the preceding [or following] layer. Being an SOS model 
means that for each column of sites along the vertical 
direction there is a unique upper occupied site, with no 
vacancies below it nor floating atoms above it. Viewed 
from above, the surface is a square network with one pair 
of diagonally opposed corners on A layers and the other 
pair on B layers. The restriction is that neighboring sites 
must be on adjacent layers (so that their separation is 
the distance from a corner to the center of the BCC lat¬ 
tice). There are then 6 possible configurations: two in 
which the two B corners are both either above or below 
the A corners and four in which one pair of catercorners 
are on the same layer and the other pair are on differ¬ 
ent layers (one above and one below the first pair). In 
the symmetric model, there are three energies, —e for the 
first pair, and ±<5/2 for the others, the sign depending on 
whether the catercorner pair on the same lattice are on 
A or B [Nolden and van Beijeren 1994]. The case <5=0 
corresponds to the E-model, which has an infinite-order 
phase transition and an essential singularity at the crit¬ 
ical point, in the class of the Kosterlitz-Thouless [1973] 
transition [Kosterlitz 1974]. (In the “ice” model, e also 
is 0.) For the asymmetric 6-vertex model, each of the 
6 configurations can have a different energy; this model 
can also be solved exactly [Yang 1967, Sutherland et al. 
1967]. 

2.2. More Formal Treatment 

To proceed more formally, we largely follow Wortis 
[1988]. The shape of a crystal is given by the length i?(h) 
of a radial vector to the crystal surface for any direction 
h. The shape of the crystal is defined as the thermo¬ 
dynamic limit of this crystal for increasing volume V, 
specifically, 

r{h,T)= \im [R{h)/aV^/^], (2) 

y—foo 

where a is an arbitrary dimensionless variable. This func¬ 
tion r(h,T) corresponds to a free energy. In particular, 
since both independent variables are field-like (and so in¬ 
trinsically intensive), this is a Gibbs-like free energy. Like 
the Gibbs free energy, r(h, T) is continuous and convex 
in h. 

The Wulff construction then amounts to a Legendre 
transformation [2] to r{h,T) from the orientation rh- 
dependent interfacial free energy fi{m, T) (or in perhaps 
more common but less explicit notation, 7 (m, T)), which 
is fp{9,T)/cos{9). For liquids, of course, /i(m,r) is 
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1 / 7 -plot for a = 0.25 


£-plot for a — 0.25 



1 / 7 -plot for a = 1.0 


^-plot for a = I.O 


FIG. 2: 7 -plots (plots of l/ 7 -plots and 7plots for 

Eq. § for positive values of a. For a = 1/4, all orientations 
appear on the ECS. For a = 1.0, the l/ 7 -plot has concave 
regions, and the has ears and flaps that must be trun¬ 

cated to give the ECS, essentially an octahedron with curved 
faces. From Sekerka [2004], which shows in a subsequent fig¬ 
ure that the 7 - and l/ 7 -plots for a — -0.2 and -0.5 resemble 
the 1 / 7 - and 7 -plots, respectively, for a = 1/4 and 1. 


spherically symmetric, as is the equilibrium shape. (Her¬ 
ring [1953] mentions rigorous proofs of this problem by 
Schwarz in 1884 and by Minkowski in 1901.) For crys¬ 
tals, /i(ih, T) is not spherically symmetric but does have 
the symmetry of the crystal lattice. For a system with 
cubic symmetry, one can write 

/7 m, T) =-fo{T)[l + a{T){m'^ + m^ + m% (3) 

where 70 (F) and a{T) are constants. As illustrated in 
Fig. i for a=l/4 the asymmetry leads to minor dis¬ 
tortions, which are rather inconsequential. However, for 
a=l the enclosed region is no longer convex, leading to 
an instability to be discussed shortly 

One considers the change in the interfacial free en¬ 
ergy associated with changes in shape. The constraint 
of constant volume is incorporated by subtracting from 
the change in the integral of /7m, T) the corresponding 
change in volume, multiplied by a Lagrange multiplier A. 



FIG. 3: Schematic of the Wulff construction. The interfacial 
free energy per nnit area bth is plotted in polar form (the 
“Wulff plot” or “ 7 -plot”). One draws a radius vector in each 
direction ih and constructs a perpendicular plane where this 
vector hits the Wulff plot. The interior envelope of the fam¬ 
ily of “Wulff planes” thus formed, expressed algebraically in 
Eq. is the crystal shape, up to an arbitrary overall scale 
factor which may be chosen as unity.” From Wortis [1988] 


Herring [1951, 1953] showed that this constrained mini¬ 
mization problem has a unique and rather simple solu¬ 
tion that is physically meaningful in the limit that it is 
satisfactory to neglect edge, corner, and kink energies in 
fi (ih, T ), that is in the limit of large volume. In this case 
A oc by choosing the proportionality constant as 

essentially the inverse of a, we can write the result as 


r(h, T) = min 

m 


/7m, T) A 

rh • h / 


(4) 


The Wulff construction is illustrated in Fig. The 
interfacial free energy at some assumed T, is dis¬ 

played as a polar plot. The crystal shape is then the 
interior envelope of the family of perpendicular planes 
(lines in 2D) passing through the ends of the radial vec¬ 
tors m/7m). Based on Eq. (|^ one can, at least in prin¬ 
ciple, determine m(h) or h(m), which thus amounts to 
the equation of state of the equilibrium crystal shape. 
One can also write the inverse of Eq. 


1 

/*(m,F) 


= min 

m 


/ i//7h,r) \ 

y ih • h j 


(5) 


Thus, a Wulff construction using the inverse of the crystal 
shape function yields the inverse free energy. 

To be more explicit, consider consider the ECS in 
Cartesian coordinates z{x,y), i.e. h oc (x,y, z(x,y)), as¬ 
suming (without dire consequences [Wortis 1988]) that 
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z{x, y) is single-valued. Then in order for any displace¬ 
ment to be tangent to h, dz — Px dx — Py dy = 0 

h = ! ^ ^{-PxZ, -Pyz, 1), (6) 

^Ji+pI+pI 

where Px is shorthand for dz/dx. 

Then the total free energy and volume are 

Fi(T) = JJ fp{px,Py)dxdy 

V = JJ z{x,y)dxdy (7) 

where fp , which incorporates the line-segment length, 
is: fp = [1 -I -+ PyV^'^ fi- Minimizing Fi subject to 
the constraint of fixed V leads to the Euler-Lagrange 
equation 

9 fp{dxZ,Py) ^ d fp(PxiPy) ro\ 

(Actually one should work with macroscopic lengths, 
then divide by the times the proportionality con¬ 

stant. Note that this leaves Px and Py unchanged. [Wor- 
tis 1988]) On the right-hand side 2A can be identified as 
the chemical potential p, so that the constancy of the 
left hand side is a reflection of equilibrium. Eq. Q is 
strictly valid only if the derivatives of fp exist, so one 
must be careful near high-symmetry orientations below 
their roughening temperature, for which facets occur. To 
show that this highly nonlinear second-order partial dif¬ 
ferential equation with unspecified boundary conditions 
is equivalent to Eq. Q, we first note that the first inte¬ 
gral of Eq. Q is simply 

z-xpx- ypy = fp{Px,Py) (9) 


H— Y —H 

• • • 



FIG. 4: Kossel crystal at T = 0. The energy to cleave the 
crystal along the depicted slanted interface (tan^ = Y/Z) is 
£1 (|A| + |y|). From Wortis [1988] 

lattice with nearest-neighbor bonds having bond energy 
ei, often called a 2D Kossel [1927, 1934] crystal; note also 
Stranski [1928]. In this model, the energy to cleave the 
crystal is the Manhattan distance between the ends of the 
cut; i.e., as illustrated in Fig.[^ the energy of severing the 
bonds between (0,0) and {X,Y) is just +ei ([A] -|- jFj). 
The interfacial area, i.e. length, is 2(A^ -|- since the 
cleavage creates two surfaces. At T = 0, entropy plays 
no role so that 

M9) = I (1 sin^l + 1 COS01) ^ I (1 + |0| + ...) (12) 


The right-hand side is just a function of derivatives, con¬ 
sistent with this being a Legendre transformation. Then 
differentiating yields 

X = -dfp/d{px), y = -dfp/d{py) ( 10 ) 

Hence, one can show that 

z{x,y) = min {fp{px, Py) + xpx + yPy) ( 11 ) 

,Py 

3. APPLICATIONS OF FORMAL RESULTS 
3.1. Cusps and Facets 

The distinguishing feature of Wulff plots of faceted 
crystals compared to liquids is the existence of [pointed] 
cusps in /i(m, T), which underpin these facets. The sim¬ 
plest way to see why the cusp arises is to examine a square 


At finite T fluctuations and attendant entropy do con¬ 
tribute, and the argument needs more care. 

Recalling Eq. 0 we see that if there is a linear cusp 
at 6 * = 0 , then 


/,(0,r) = /,(o,r) + H(T)i0i, 


(13) 


where B = (3{T)/h, since the difference between / (9) 
and fp{9) only appears at order 9“^. Comparing Eqs. (12) 
and (13), we see that for the Kossel square fi{0, 0) = eij2 
and BJO) = ei/2. Further discussion of the 2D fi{9) is 
deferred to Section 4.3 below. 

To see how a cusp in fi{m,T) leads to a facet in the 
ECS, consider Fig.l^ the Wulff plane for 0 > 0 intersects 
the horizontal 9 = O^plane at a distance fi{0) + d{9) from 
the vertical axis. The crystal will have a horizontal axis 
if and only if d{9) does not vanish as 0 —0. From Fig. 
it is clear that 9 « sin0 « B9/d{9) for 9 near 0, so that 
d{0) = B > 0. For a weaker dependence on 9, e.g. B\9\‘' 
with C > 1, d(0) = 0, and there is no facet. Likewise, 
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FIG. 5: Wulff plot with a linear cusp at 0 = 0. If d{6) —>■ 0 
as 0 —>■ 0, then there is no facet corresponding to 0 = 0, and 
the 0 = 0 Wulff plane (dotted line) is tangent to the crystal 
shape at just a single point. Since d{6) = B, a cusp in the 
Wulff plot leads to a facet of the corresponding orientation 
on the ECS. From Wortis [1988] 


at the roughening temperature ^ vanishes and the facet 
disappears. 


3.2. Sharp Edges and Forbidden Regions 

When there is a sharp edge (or corner) on the ECS 
r(h, T), Wulff planes with a range of orientations rh will 
not be part of the inner envelope determining this ECS; 
they will lie completely outside it. There is no portion of 
the ECS whose surface tangent has these orientations. As 
in the analogous problems with forbidden values of the 
“density” variable, the free energy /i(m,T) is actually 
not properly defined for forbidden values of m; those 
unphysical values should actually be removed from the 
Wulff plot. Fig. [^depicts several possible ECSs and their 
associated Wulff plots. It is worth emphasizing that, in 
the extreme case of the fully faceted ECS at T = 0, the 
Wulff plot is simply a set of discrete points in the facet 
directions. 

Now if we denote by m+ and m_ the limiting orien¬ 
tations of the tangent planes approaching the edge from 
either side, then all intermediate values do not occur as 
stable orientations. These missing, not stable, “forbid¬ 
den” orientations are just like the forbidden densities at 
liquid-gas transitions, forbidden magnetizations in ferro- 
magnets at T < [Garcia 1984], and miscibility gaps 
in binary alloys. Herring [1951, 1953] first presented an 
elegant way to determine these missing orientations us¬ 
ing a spherical construction. For any orientation rh, this 
tangent sphere (often called a “Herring sphere”) passes 
through the origin and is tangent to the Wulff plot at 
/i(ih). From geometry he invoked the theorem that an 
angle inscribed in a semicircle is a right angle. Thence, 
if the orientation rh appears on the ECS, it appears at 



FIG. 6: Some possible Wulff plots and corresponding equilib¬ 
rium crystal shapes. Faceted and curved surfaces may appear, 
joined at sharp or smooth edges in a variety of combinations. 
From Rottman and Wortis [1984a]; the EGS are also in Her¬ 
ring [1953] 


an orientation that points outward along the radius of 
that sphere. Herring then observes that only if such a 
sphere lies completely inside the plot of /i(m) does, that 
orientation appear on the ECS. If some part were inside, 
its Wulff plane would clip off the orientation of the point 
of tangency, so that that orientation would be forbidden. 

The origin of a hill-and-valley structure from the con¬ 
stituent free energies [Jeong and Williams 1999, Williams 
and Barteltl996] is illustrated schematically in Fig. It 
arises when they satisfy the inequality 

/j(m = ni)Ai -b fi{n 2 )A 2 < /i(n)A, (14) 

where Ai and A 2 are the areas of strips of orientation 
Ill and n 2 , respectively, while A is the area of the sum 
of these areas projected onto the plane bounded by the 
dashed lines in the figure. This behavior, again, is con¬ 
sistent with the identification of the misorientation as 
a density (or magnetization)-like variable rather than a 
field-like one. 
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It can then be shown that m[/i(m)]avr lies on the de¬ 
picted circle, so that the Wulff plane passes through the 
edge at R. 


3.3. Going Beyond Wulff Plots 


FIG. 7: Illustration of how orientational phase separation oc¬ 
curs when a “hill-and-valley” structure has a lower total sur¬ 
face free energy per area than a flat surface as in Eq. (14l. 
The sketch of the free energy vs. p = taxi 6 shows that this 
situation reflects a region with negative convexity which is 
accordingly not stable. The dashed line is the tie bar of a 
Maxwell or double-tangent construction. The misorientations 
are the coexisting slat-like planes, with orientations ni and 
02 , in the hill-and-valley structure. From Jeong and Williams 
[1999] 


The details of the lever rule for coexistence regimes 
were elucidated by Wortis [1988]: As depicted in Fig. 
which denotes as P and Q the two orientations bounding 
the region that is not stable, the lever-rule interpolations 
lie on segments of a spherical surface. Let the edge on 
the ECS be at R. Then an interfaced created at a forbid¬ 
den rh will evolve toward a hill-and-valley structure with 
orientations P and Q with a free energy per area of 


[/i(lSl)]avr 


x/.(P) + y/»(Q) 

d 


(15) 



FIG. 8: ECS analogue of the Maxwell double-tangent con¬ 
struction. O is the center of the crystal. Points P and Q 
are on the (stable) Wulff plot but the region between them 
is unstable; hence, the ECS follows PRQ and has an edge at 
R. An interface at the intermediate orientation ih breaks up 
into the orientations P and Q with relative proportions x:y ; 
thus, the average free energy per unit area given by Eq. (151, 
which in turn shows that /i(m)]avr lies on the circle. From 
Wortis [1988] 


To determine the limits of forbidden regions, it is more 
direct and straightforward to carry out a polar plot of 
l//i(m) [Frank, 1963] rather than /i(m), as discussed 
in Sekerka’s [2004] review chapter. Then a sphere pass¬ 
ing through the origin becomes a corresponding plane; 
in particular, a Herring sphere for some point becomes 
a plane tangent to the plot of l//i(m). If the Herring 
sphere is inside the Wulff plot, then its associated plane 
lies outside the plot of l//i(m). If, on the other hand, 
if some part of the Wulff plot is inside a Herring sphere, 
the corresponding part of the l//i (ih) plot will be outside 
the plane. Thus, if the plot of l//i(m) is convex, all its 
tangent planes will lie outside, and all orientations will 
appear on the ECS. If it is not convex, it can be made 
so by adding tangent planes. The orientations associated 
with such tangent planes are forbidden, so their contact 
curve with the l//i(m) plot gives the bounding stable 
orientations into which forbidden orientations phase sep¬ 
arate. 

Summarizing the discussion in Sekerka [2004], the con¬ 
vexity of l//i(m) can indeed be determined analyti¬ 
cally since the curvature l//i(m) is proportional (with 
a positive-definite proportionality constant) to the stiff¬ 
ness, i.e. in 2D, 7 -|- = 7 , or preferably (3 + 

d^P/dO'^ = ,5 as in Eq. (Ilj) to emphasize that the stiff¬ 
ness and [step] free energy per length have different units 
in 2D from 3D. Hence, l//i(rh) is not convex where the 
stiffness is negative. The very complicated generaliza¬ 
tion of this criterion to 3D is made tractable via the 
^-vector formalism of Hoffman and Cahn [1972, 1974], 
where ^ = V(r/i(m)), where r is the distance from the 
origin of the 7 plot. Thus, 


/j(m)=|-m, m-d| = 0 , (16) 

which is discussed well by Wheeler [1999] and Sekerka 
[2004]. To elucidate the process, we consider just the 2D 
case [Cahn and Carter 1996]; cf. Fig.[^ 

The solid curve in Fig. is the | plot and the dashed 
curve is the l/ 7 -plot for /i(m) = 7 oc 1 -|- 0.2 cos 40. For 
this case the 1 / 7 -plot is not convex and the ^ plot forms 
“ears.” The equilibrium shape is given by the interior 
envelope of the ^ plot; in this case it exhibits four corners. 

Pursuing this analogy, we see that if one cleaves a crys¬ 
tal at some orientation rh that is not on the ECS, i.e. 
between th^ and rh_, then this orientation will break up 
into segments with orientations rh+ and rh_ such that 
the net orientation is still rh, providing another example 
of the lever rule associated with Maxwell double-tangent 
constructions for the analogous problems. The time to 












FIG. 9: Graphical constructions for an anisotropic fi{m) for 
various values of an anisotropy parameter a, where /; oc 
1 + a cos^ 6 ^ sin^ In the left column fi{9) is plotted from 
top to bottom for a = 1/2, 1, 2. Anisotropy increases with 
positive a, so 1 /a corresponds in some sense to a temperature 
in conventional plots. In the center panel, rii is cos^ 6. The 
shape resulting from the gradient construction with the ears 
removed is the Wulff EGS. From Cahn and Carter [1996] 



FIG. 10: The solid curve is the ^ plot, while the dashed curve 
is the 1 / 7 -plot for fi{m) = 7 oc 1 -I- 0.2 cos40. For this case 
(but not for small values of a) the l/ 7 -plot is not convex, and 
the ^ plot forms “ears.” These ears are then removed, so that 
the equilibrium shape is given by the interior envelope of the 
^ plot, in this case having four corners. From Wheeler [1999] 


evolve to this equilibrium state depends strongly on the 
size of energy barriers to mass transport in the crystalline 
material; it could be exceedingly long. To achieve rapid 
equilibration, many nice experiments were performed on 
solid hep ^He bathed in superfluid ^He, for which equili¬ 
bration occurs in seconds (Balibar and Castaing [1980], 
Keshishev et al. [1981], Wolf et al. [1983, 1985], and many 
more; see Balibar et al. [2005] for a comprehensive recent 
review. Longer but manageable equilibration times are 
found for Si and for An, Pb, and other soft transition 
metals. 


4. SOME PHYSICAL IMPLICATIONS OF 
WULFF CONSTRUCTIONS 

4.1. Thermal Faceting and Reconstruction 

A particularly dramatic example is the case of surfaces 
vicinal to Si (111) by a few degrees. In one misorien- 
tation direction the vicinal surface is stable above the 
reconstruction temperature of the (111) facet but below 
that temperature, /^(lll) decreases significantly so that 
the original orientation is no longer stable and phase sep¬ 
arates into reconstructed (111) terraces and more highly 
misoriented segments [Phaneuf, 1987; Bartelt, 1989]. 
The correspondence to other systems with phase separa¬ 
tion at first-order transitions is even more robust. Within 
the coexistence regime one can in mean field determine a 
spinodal curve. Between it and the coexistence boundary 
one observes phase separation by nucleation and growth, 
as for metastable systems; inside the spinodal one ob¬ 
serves much more rapid separation with a characteristic 
most-unstable length [Phaneuf, 1993]. This system is dis¬ 
cussed further in Section 7.1 below. Furthermore, there 
are remarkably many ordered phases at larger misori- 
entations [Olshanetsky and Mashanov 1981, Baski and 
Whitman 1997]. 

Wortis [1988] describes “thermal faceting” experiments 
in which metal crystals, typically late transition or noble 
metal elements like Cu, Ag, and Fe, are cut at a high 
Miller-index direction and polished. They are then an¬ 
nealed at high temperatures. If the initial plane is in a 
forbidden direction, optical striations, due to hill and val¬ 
ley formation, appear once these structures have reached 
optical wavelengths. While the characteristic size of this 
pattern continues to grow as in spinodal decomposition, 
the coarsening process is eventually slowed and halted by 
kinetic limitations. 

There are more recent examples of such phenomena. 
After sputtering and annealing above 800K, Au(4,5,5) at 
300K forms a hill-and-valley structure of two Au(lll) 
vicinal surfaces, one that is reconstructed and the other 
which is not, as seen in Fig. El This seems to be an equi¬ 
librium phenomenon: it is reversible and independent of 
cooling rate [Rousset 1999]. Furthermore, while it has 
been long known that adsorbed gases can induce faceting 
on bcc (111) metals [Bonezek et al. 1980], ultrathin metal 
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FIG. 11: Morphology of the faceted Au(4,5,5) surface mea¬ 
sured at room temperature, (a) 3D plot of a large-scale (scan 
area: 1.4 mm x 1.4 mm) STM image. Phases A and B form 
the hill-and-valley morphology, (b) STM image zoomed in on 
a boundary between the two phases. All steps single-height, 
i.e. 2 . 35 A high. Phase B has smaller terraces, 13A wide, while 
phase A terraces are about 30A wide. This particular surface 
has (2,3,3) orientation. From Rousset [1999] 


films have also been found to produce faceting of W(lll), 
W(211), and Mo(lll) [Madey et al. 1996, 1999]. 


4.2. Types A &: B 

The above analysis indicates that at T = 0 the ECS 
of a crystal is a polyhedron having the point symme¬ 
try of the crystal lattice, a result believed to be general 
for finite-range interactions [Fisher 1983]. All boundaries 
between facets are sharp edges, with associated forbid¬ 
den non-facet orientation; indeed, the Wulff plot is just 
a set of discrete points in the symmetry directions. At 
finite temperature, two possibilities have been delineated 
(with cautions [Wortis 1988], labeled (nonmnemonically) 
A and B. In type-A, there are smooth curves between 
facet planes rather than edges and corners. Smooth here 
means, of course, that not only is the ECS continuous, 
but so is its slope, so that there are no forbidden orienta¬ 
tions anywhere. This situation corresponds to continuous 
phase transitions. In type B, in contrast, corners round 
at finite T but edges stay sharp until some temperature 



FIG. 12: Interfacial phase diagrams for simple-cubic nearest- 
neighbor Kossel crystal with nearest-neighbor as well as 
(weak) next-nearest-neighbor (NNN) attractions. The angu- 
lar variables 9 and <j) (not to be confused with (p, cf. Section 


2 2.11 interfacial orientation (m) and EGS (h), respectively. 


in an equatorial section of the full 3D phase diagram, a) 
The T — 6 phase diagram (a) shows the locus of cusps in the 
Wulff plot along the symmetry directions below the respective 
roughening temperatures. For no NNN interaction (e 2 = 0), 
there are only cusps at vertical lines at 0 and 7r/2. b) The 
T — h phase diagram gives the faceted areas of the crystal 
shape. The NNN attraction leads to additional (111) [not 
seen in the equitorial plane] and (110) facets at low enough 
temperature. Thus, for £2 = 0 the two bases of the (100) 
and (010) phases meet and touch each other at (and only 
at) 0 = 7r/4 (at T = 0), with no intervening (110) phase. 
Each type of facet disappears at its own roughening tempera¬ 
ture. Above the phase boundaries enclosing those regions, the 
crystal surfaces are smoothly curved (i.e., thermodynamically 
“rough”). This behavior is consistent with the observed phase 
diagram of hep ‘*He. From Rottman and Wortis ]1984a,b] 


Tq. For To < T < Ti there are some rounded edges and 
some sharp edges, while above Ti all edges are rounded. 

Rottman and Wortis [1984a,b] presents a comprehen¬ 
sive catalogue of the orientation phase diagrams, Wulff 
plots, and ECSs for the cases of non-existent, weakly 
attractive, and weakly repulsive next-nearest-neighbor 
(NNN) bonds in 3D. Figs. [T^ and 13 show the orienta¬ 
tion phase diagrams and the Wulff plots with associated 
ECSs, respectively, for weakly attractive NNN bonds. As 
indicated in the caption, it is easy to describe what then 
happens when £2=0 and only {100} facets occur. Like¬ 
wise, Figs.fT^and 15 show the orientation phase diagrams 
and the Wulff plots with associated ECSs, respectively, 
for weakly repulsive NNN bonds. 


4.3. 2D Studies 


Exploring the details is far more transparent in 2D 
than 3D. The 2D case is physically relevant in that it 
describes the shape of islands of atoms of some species 
at low fractional coverage on an extended flat surface of 
the same or another material. An entire book is devoted 
to 2D crystals [Lyuksyutov et al 1992]. The 2D perspec¬ 
tive can also be applied to cylindrical surfaces in 3D, as 
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Wulff 


T-0 


0<T<T, 



Crystal Shape 
(equatorial (perspective) 

plane) 






FIG. 13: Representative Wulff plots and ECS’s for the crystal 
with weak NNN attractions whose phase diagram is shown in 
Fig |12| At low enough temperature there are (100), (110) and 
(111) facets. For weak attraction the (110) and (111) facets 
roughen away below the (100) roughening temperature. For 
£2 = 0, Ti {2 = 0, so that the configurations in the second row 
do not occur; in the first row, the octagon becomes a square 
and the perspective shape is a cube. Facets are separated at 
T > 0 by curved surfaces, and all transitions are second order. 
Spherical symmetry obtains as T approaches melting at Tc. 
From Rottman and Wortis [1984] 


remarked by Nozieres [1992]. Formal proof is also more 
feasible, if still arduous, in 2D; an entire book is devoted 
to this task [Dobrushin et al. 1992]; see also Pfister [1991] 
and Miracle-Sole and Ruiz [1994]. 

For the 2D nearest-neighbor Kossel crystal described 
above Wortis [1988] notes that at T = 0 a whole class of 
Wulff planes pass through a corner. At finite T thermal 
fluctuations lift this degeneracy and the corner rounds, 
leading to type A behavior. To gain further insight, we 
now include a next-nearest-neighbor (NNN) interaction 
£ 2 , so that 

f^i^) = (|cos6l| + |sin6»|) -b ^ (|cos6l| - |sin6»|) 

(17) 

For favorable NNN bonds, i.e. €2 > 0, one finds new {11} 
facets but still type-A behavior with sharp edges, while 
for unfavorable NNN bonds, i.e. £2 < 0, there are no new 
facets but for finite T the edges are no longer degenerate 
so that type-B behavior obtains. Again recalling that 
fi{9) = fp{6) I cos0|, we can identify /g = £2 -b ei/2 and 
P/h = £ 1 / 2 , as noted in other treatments, e.g. Dieluweit 
et al. [2003]. That work, however, finds that such a 



FIG. 14: Interfacial phase diagram with (weak) next-nearest- 


neighbor (NNN) repulsion rather than attraction as in Fig. 12 


The NNN repulsion stabilizes the (100) facets. Gurved sur¬ 
faces first appear at the cube corners and then reach the equa¬ 
torial plane at Ta. The transition at the equator remains first 
order until a higher temperature Tt. The dotted boundaries 
are first order. A forbidden (coexistence) region appears in 
the T—h phase diagram. From Rottman and Wortis [1984a,b] 


Wulff 

Plot 


0£T<T3 


Crystal Shope 

(equatorial (perspective) 
plane) 



T3<T<T, 








Tr<T<T, 




FIG. 15: Representative Wulff plots and EGS’s for the crys¬ 
tal with weak NNN repulsions whose phase diagram is shown 
in Fig[g Gurved surfaces appear hrst at the cube corners. 
Junctions between facets and curved surfaces may be either 
first- or second-order (sharp or smooth), depending on orien¬ 
tation and temperature. From Rottman and Wortis [1984a,b] 


model cannot adequately account for the orientation- 
dependent stiffness of islands on Cu(OOl). Attempts to 
resolve this quandry using 3-site non-pairwise (trio) in¬ 
teractions [Stasevich et al. 2004, 2006] did not prove en¬ 
tirely satisfactory. In contrast, on the hexagonal Cu(lll) 
surface, only NN interactions are needed to account ade¬ 
quately for the experimental data [Stasevich et al. 2005, 





































































11 


2006]. In fact, for the NN model on a hexagonal grid, Zia 
[1986] found an exact and simple, albeit implicit, exact 
expression for the ECS. However, on such (111) surfaces 
(and basal planes of hep crystals) lateral pair interac¬ 
tions alone cannot break the symmetry to produce a dif¬ 
ference in energies between the two kinds of step edges, 
viz. {100} and {111}) microfacets (A and B steps, re¬ 
spectively, with no relation to types A and B!). The 
simplest viable explanation is an orientation-dependent 
trio interaction; calculations of such energies support this 
idea [Stasevich et al. 2005, 2006]. 

Strictly speaking, of course, there should be no 2D 
facet (straight edge) and accompanying sharp edges (cor¬ 
ners) at r > 0 [Gallavotti 1972; Abraham and Reed 1974, 
1976; Avron et al. 1982 and references therein] since that 
would imply ID long-range order, which should not oc¬ 
cur for short-range interactions. Measurements of islands 
at low temperatures show edges that appear to be facets 
and satisfy Wulff corollaries such as that the ratio of the 
distances of two unlike facets from the center equals the 
ratio of their fi [Kodambaka et al. 2006]. Thus, this is¬ 
sue is often just mentioned in passing [Michely and Krug, 
2004] or even ignored. On the other hand, sophisticated 
approximations for fi(9) for the 2D Ising model, includ¬ 
ing NNN bonds have been developed, e.g. [Stasevich 
and Einstein 2007], allow numerical tests of the degree 
to which the ECS deviates from a polygon near corners 
of the latter. One can also gauge the length scale at which 
deviations from a straight edge come into play by using 
that the probability per atom along the edge for a kink 
to occur is essentially the Boltzmann factor associated 
with the energy to create the kink [Weeks 2014]. 

Especially for heteroepitaxial island systems (when the 
island consists of a different species from the substrate), 
strain plays an important if not dominant role. Such sys¬ 
tems have been investigated, e.g., by E. Liu [2006], who 
points out that for such systems the shape does not sim¬ 
ply scale with A, presumably implying the involvement 
of some new length scale [s]. A dramatic manifestation 
of strain effects is the island shape transition of Cu on 
Ni(OOl), which changes from compact to ramified as is¬ 
land size increases [Muller et al. 1998]. Eor small islands, 
additional quantum-size and other effects lead to favored 
island sizes (magic numbers). 


5. VICINAL SURFACES-ENTREE TO ROUGH 
REGIONS NEAR FACETS 


In the rough regions the ECS is a vicinal surface of 
gradually evolving orientation. To the extent that a 
local region has a particular orientation, it can be ap¬ 
proximated as an infinite vicinal surface. The direction 
perpendicular to the terraces (which are densely-packed 
facets) is typically called z. In “Maryland notation” (cf. 
[|2 2.1 )the normal to the vicinal surface lies in the x — z 
plane, and the distance ^ between steps is measured along 
i, while the steps run along the y direction. In the sim¬ 


plest and usual approximation, the repulsions between 
adjacent steps arise from two sources: an entropic or 
steric interaction due to the physical condition that the 
steps cannot cross, since overhangs cannot occur in na¬ 
ture. The second comes from elastic dipole moments due 
to local atomic relaxation around each step, leading to 
frustrated lateral relaxation of atoms on the terrace plane 
between two steps. Both interactions are oc 1/P. 

The details of the distribution P{t) of spacings between 
steps have been reviewed in many places [Jeong and 
Williams 1999, Einstein et al. 2001, Giesen 2001, Einstein 
2007] The average step separation {tj is the only charac¬ 
teristic length in the x direction. N.B., {t) need not be a 
multiple of, or even simply related to, the substrate lat¬ 
tice spacing. Therefore, we consider P(s) = (£)“^P(£), 
where s = ^/(^), a dimensionless length. Eor a “per¬ 
fect” cleaved crystal, P[s) is just a spike (5(s — 1). For 
straight steps placed randomly at any position with prob¬ 
ability 1/(1'), P{s) is a Poisson distribution exp(—s). Ac¬ 
tual steps do meander, as one can study most simply 
in a terrace-step-kink (TSK) model. In this model, the 
only excitations are kinks (with energy e) along the step. 
(This is a good approximation at low temperature T since 
adatoms or vacancies on the terrace cost several e\ [4ei 
in the case of a simple cubic lattice].) The entropic re¬ 
pulsion due to step meandering dramatically decreases 
the probability of finding adjacent steps at 1 {€)■ To 

preserve the mean of one, P{s) must also be smaller than 
exp(—s) for large s. 

If there is an additional energetic repulsion A/P, the 
magnitude of the step meandering will decrease, nar¬ 
rowing P{s). As A —>■ 00 , the width approaches 0 
(P(s) —> 5{s — 1), the result for perfect crystals). We 
emphasize that the energetic and entropic interactions 
do not simply add. In particular, there is no negative 
(attractive) value of A at which the two cancel each 
other. (Cf. Eq. (27) below.) Thus, for strong repulsions, 
steps rarely come close, so the entropic interaction plays 
a smaller role, while for A < 0, the entropic contribution 
increases, as illustrated in Fig. [T6|and explicated below. 
We emphasize that the potentials of both interactions 
decay as £~^ (cf. Eq. (27 below), in contrast to some 
statements in the literature (in papers analyzing ECS 
exponents) that entropic interactions are strictly short 
range while energetic ones are long-range. 

Investigation of the interaction between steps has been 
reviewed well in several places [Jeong and Williams 1999, 
Giesen 2001, Nozieres 1999, Einstein 1996, Einstein et 
al. 2001]. The earliest studies seeking to extract A from 
terrace-width distributions (TWDs) used the mean-field¬ 
like Gruber-Mullins [Gruber and Mullins 1967] approxi¬ 
mation, in which a single active step fluctuates between 
two fixed straight steps 2(£} apart. Then the energy as¬ 
sociating with the fluctuations x{y,t) is 


Af = -mLy + dy, (18) 

where Ly is the size of the system along the mean step 
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direction (i.e. the step length with no kinks). We ex¬ 
pand /3(d) as the Taylor series /3(0) -I- /3'(0)d -I-1/2 /3"(0)d^ 
and recognize that the length of the line segment has 
increased from dy to dy /cos 6 Ri dy{l -I- 1/2 0^). For close- 
packed steps, for which /3'(0) = 0, it is well known that 
(using 9 « tan 6 = dx/dy) 


AS 


m 

2 


dy 


dy, /3(0) = /3(0)+/3"(0), (19) 


where /3 is the step stiffness [Fisher et al. 1982]. N.B., 
the stiffness /3(d) has the same definition for steps with 
arbitrary in-plane orientation — for which /3'(d) 7 ^ 0 — 
because to create such steps, one must apply a “torque” 
[Leamy et al. 1975] which exactly cancels /3'(d). (See 
Stasevich [2006, 2007] for a more formal proof.) 

Since x{y) is taken to be a single-valued function that 
is defined over the whole domain of y, the 2D configu¬ 
ration of the step can be viewed as the worldline of a 
particle in ID by recognizing y as a time-like variable. 
Since the steps cannot cross, these particles can be de¬ 
scribed as spinless fermions in ID, as pointed out first 
by de Gennes [1968] in a study of polymers in 2D. Thus, 
this problem can be mapped into the Schrodinger equa¬ 
tion in ID: dx/dy in Eq. (19) becomes dx/dt, with the 


form of a velocity, with the stiffness playing the role of 
an inertial mass. This correspondence also applies to 
domain walls of adatoms on densely-covered crystal sur¬ 
faces, since these walls have many of the same properties 
as steps. Indeed, there is a close correspondence between 
the phase transition at smooth edges of the ECS and the 
commensurate-incommensurate phase transitions of such 
overlayer systems, with the rough region of the ECS cor¬ 
responding to the incommensurate regions and the local 
slope related to the incommensurability [Pokrovsky and 
Talapov 1979, 1984, Villain 1980, Haldane and Villain 
1981, Schulz et al. 1982]. Jayaprakash et al. [1984] pro¬ 
vide the details of the mapping from a TSK model to 
the fermion picture, complete with fermion creation and 
annihilation operators. 

In the Gruber-Mullins [1967] approximation, a step 
with no energetic interactions becomes a particle in a 
ID infinite-barrier well of width 2{t}, with well-known 
ground-state properties 


{TTkeTf 


oc sin f 77 ^^ ; P(s) = sin^ (; Eg = 

^ { 2 J’ 8/3(^)2 

. 

Thus, it is the kinetic energy of the ground state in the 
quantum model that corresponds to the entropic repul¬ 
sion (per length) of the step. In the exact solution for 
the free energy expansion of the equilibrium crystal shape 
[Akutsu et al. 1988], the numerical coefficient in the cor¬ 
responding term is i/e rather than i/s. Note that P(s) 
peaks at s = 1 and vanishes for s > 2 . 

Suppose, next, that there is an energetic repulsion 
[/(£) = A/S^ between steps. In the ID Schrodinger equa¬ 


tion, the prefactor of —d‘^i\)(S)ldf' is (/c_bF)^/ 2,5, with 
the thermal energy ksT replacing h. (Like the repulsions, 
this term has units £~'^.) Hence, A only enters the prob¬ 
lem in the dimensionless combination A = Ad/iksT)'^ 
[Jeong and Weeks 1999]. In the Gruber-Mullins picture, 
the potential (per length) experienced by the single ac¬ 
tive particle is (with £ = £—{£)): 


U{£) = 


A 


A 


(£-(£}y (£+(£}r 



The first term is just a constant shift in the energy. For 
A sufficiently large, the particle is confined to a region 
\£\ <C {£), so that we can neglect the fixed walls and the 
quartic term, reducing the problem to the familiar simple 
harmonic oscillator, with the solution: 


V'o(^)oc e = 






exp 


2(Tq 


( 22 ) 

where ac = (48A) and w = 

For A of 0 or 2, the TWD can be computed exactly 
(See below). For these cases, Eqs. (20) and (22), respec¬ 
tively, provide serviceable approximations. It is Eq. (22) 


that is prescribed for analyzing TWDs in the most-cited 
resource on vicinal surfaces [Jeong and Williams 1999]. 
Indeed, it formed the basis of initial successful analyses 
of experimental STM (scanning tunneling microscopy) 
data [Wang et al. 1990]. However, it has some notable 
shortcomings. Perhaps most obviously, it is useless for 
small but not vanishing A, for which the TWD is highly 
skewed, not resembling a Gaussian, and the peak, corre¬ 
spondingly, is significantly below the mean spacing. For 
large values of A, it significantly underestimates the vari¬ 
ance or, equivalently, the value of A one extracts from 
the experimental TWD width [Ihle et al.[1998]]: in the 
Gruber-Mullins approximation the TWD variance is the 
same as that of the active step, since the neighboring step 
is straight. For large A the fluctuations of the individ¬ 
ual steps on an actual vicinal surface become relatively 
independent, so the variance of the TWD is the sum of 
the variance of each, i.e. twice the step variance. Given 
the great (quartic) sensitivity of A to the TWD width, 
this is problematic. As experimentalists acquired more 
high-quality TWD data, other approximation schemes 
were proposed, all producing Gaussian distributions with 
widths (X but with proportionality constants no¬ 

tably larger than 48“^/"* = 0.38. 

For the “free-fermion” {A = 0) case. Jobs et al. [1991] 
developed a sequence of analytic approximants to the ex¬ 
act but formidable expression [Dyson 1962, Mehta 2004] 
for P{s). They as well as a slightly earlier paper [Bartelt 
et al. 1990] draw the analogy between the TWD of vic¬ 
inal surfaces and the distribution of spacings between 
interacting (spinless) fermions on a ring, the Calogero- 
Sutherland model Calogero 1969, Sutherland 1971, which 
in turn for three particular values of the interaction— 
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in one case repulsive {A = 2), in another attractive 
{A = — V 4 ), and lastly the free-fermion case {A = 0)— 
could be solved exactly by connecting to random matrix 
theory [Mehta 2004, Dyson 1970, Einstein 2007]; Fig. 5 
of [Jobs 1991] depicts the three resulting TWDs. 

These three cases can be well described by the Wigner 
surmise, for which there are many excellent reviews [Guhr 
et al. 1998, Mehta 2004, Haake 1991]. Explicitly, for 
p = 1, 2, and 4, 


Pg{s) = OgS^exp {-hgS^) , (23) 

where the subscript of P refers to the exponent of s. In 
random matrix literature, the exponent of s, viz. 1, 2, or 
4, is called j5, due to an analogy with inverse temperature 
in one justification. However, to avoid possible confusion 
with the step free energy per length /3 or the stiffness /3 
for vicinal surfaces, I have sometimes called it instead by 
the Greek symbol that looked most similar, g, and do so 
in this chapter. The constants bg, which fixes its mean 
at unity, and Ug, which normalizes P{s), are 


- r(^) 

(24) 

Specifically, bg = tt/4, 4 / 7 r, and 64/97r, respectively, 
while Ug = 1112, ?>2 Iit'^, and (64/97r)^, respectively. 

As seen most clearly by explicit plots, e.g. Fig. 4.2a of 
Haake [1991], T’i(s), ^ 2 ( 3 ), and Pi{s) are excellent ap¬ 
proximations of the exact results for orthogonal, unitary, 
and symplectic ensembles, respectively, and these simple 
expressions are routinely used when confronting experi¬ 
mental data in a broad range of physical problems [Guhr 
et al. 1998, Haake 1991]. (The agreement is particularly 
outstanding for P 2 {s) and ^ 4 ( 5 ), which are the germane 
cases for vicinal surfaces, significantly better than any 
other approximation [Gebremariam 2004]. 

Thus, the Calogero-Sutherland model provides a con¬ 
nection between random matrix theory, notably the 
Wigner surmise, and the distribution of spacings between 
fermions in ID interacting with dimensionless strength A. 
Specifically, 

i=f(f-l) ^ g=l + Vl + 4A. (25) 



For an arbitrary system, there is no reason that A should 
take on one of the three special values. Therefore, we 
have used Eq. (251 for arbitrary g or A, even though there 
is no symmetry-based justification of distribution based 
on the Wigner surmise of Eq. (23), and refer hereafter 
to this formula Eqs. ( 23|24 ) as the GWD (generalized 
Wigner distribution). Arguably the most convincing ar¬ 
gument is a comparison of the predicted variance with 
numerical data generated from Monte Carlo simulations. 
See Einstein [2007] for further discussion. 

There are several alternative approximations that lead 
to a description of the TWD as a Gaussian. Ihle, Pierre- 
Louis, and Misbah [1998], in particular, focus on the 
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FIG. 16: Illustration of how entropic repulsion and energetic 
interactions combine, plotted vs. the dimensionless energetic 
interaction strength A = AP/{kBT)^. The dashed straight 
line is just A. The solid curve above it is the combined en¬ 
tropic and energetic interactions, labeled Aeff for reasons ex¬ 
plained below. The difference between the two cnrves at any 
valne of the abscissa is the dimensionless entropic repulsion for 
that A. The decreasing cnrve, scaled on the right ordinate, is 
the ratio of this entropic repulsion to the total dimensionless 
repulsion AeS. It falls monotonically with A, passing thro ugh 
unity at A = 0. See the discussion accompanying Eq. (26) 
for more information and explicit expressions for the curves. 
From Ei07. 


limit of large g, neglecting the entropic interaction in 
that limit. The variance oc A~^l'^, the proportional¬ 
ity constant is 1.8 times that in the Gruber-Mullins case. 
This approximation is improved, especially for repulsions 
that are not extremely strong, by including the entropic 
interaction in an average way. This is done by replacing 
iby 


i.,= (|)Gi + |. (26) 

Physically, Aefi gives the full strength of the inverse- 
square repulsion between steps, i.e. the modification due 
to the inclusion of entropic interactions. Thus, in Eq. ([^ 


9{T) 


{nksT)'^ 


AcS 


{irksTY 

24/i3/3 


1 -b \/l+4i 


■ (27) 


From Eq. (26) it is obvious that the contribution of 
the entropic interaction, viz. the difference between the 
total and the energetic interaction, as discussed in con¬ 
junction with Fig. is p/2. Remarkably, the ratio 
of the entropic interaction to the total interaction is 
(p/2)/(p/2)^ = 2/p; this is the fractional contribution 
that is plotted in Fig. [T^ 
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FIG. 17: Critical behavior of the crystal shape near a smooth 
(second-order) edge, represented by the dot at (xo,zo)- The 
temperature is lower than the roughening temperature of the 
facet orientation, so that the region to the left of the dot is 
flat. The curved region to the right of the dot correspond to 
a broad range of rough orientations. In the thermodynamic 
limit, the shape of the smoothly curved region near the edge 
is described by the power law z ^ zq — {x — xo)^■ Away from 
the edge there are “corrections to scaling”, i.e. higher-order 
terms (cf. Eq. (331. For an actual crystal of any finite size, 
there is “finite-size rounding” near the edge, which smooths 
the singular behavior. Adapted from Jayaprakash and Saam 
[1984] 


6. CRITICAL BEHAVIOR OF ROUGH 
REGIONS NEAR FAGETS 

6.1. Theory 


Assuming (cf. Fig. 17) z the direction normal to the 
facet and {xq, zq) denote the facet edge, z ~ Zo — {x — Xo)'^ 
for X > xq. We show that the critical exponent m has 
the value 3/2 for the generic smooth edge described by 
Eq. 0 (with the notation of Eq. 


fp{p) = fo +Bp + gp^ + cp^. (28) 

Then we perform a Legendre transformation [Callen, 
1985] as in Andreev [1981] (foreshadowed in Landau and 
Lifshitz 1980) and Jayaprakash and Saam [1984]; explic¬ 
itly 


fpip) - fiv) 


dSj 

dp 


^ = rj = B + 3gp^ + 4cp^ (29) 


Hence 


fiv) = fo- ‘^gp^iv) - 3c/(ry) (30) 


But from Eq. (29) 




^ 2c f g—B 


35 V 3g 


1/2 


(31) 


Inserting this into Eq. (30) gives 


(32) 

ioi g > B and f{g) = fo for g < B. (See Jayaprakash 
et al. [1983]; Jayaprakash and Saam [1984], van Beijeren 
and Nolden [1987].) Note that this result is true not just 
for the free fermion case but even when steps interact. 
Jayaprakash et al. [1984] further show that the same d 
obtains when the step-step interaction decreases with a 
power law in i that is greater than 2. We identify f{g) 
with r(h), i.e. the magnetic-field-like variable discussed 
corresponds to the so-called Andreev field g. Writing 
Zo = /o/^ and xq = B/X, we find the shape profile 


z{x) 

zo 


= 1-2 


„ \ 1/2 / X 3/2 

Jo\ I X-Xo\ 


Cfo / X-Xq 

5^ V ^0 
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-) 


X — Xo 
Zo 


5/2 


(33) 


Note that the edge position depends only on the step free 
energy B, not on the step repulsion strength; conversely, 
the coefficient of the leading {x — xo)^''^ term is indepen¬ 
dent of the step free energy but varies as the inverse root 
of the total step repulsion strength, i.e. as g~^l'^■ 

If instead of Eq. (28) one adopts the phenomenologi¬ 


cal Landau theory of continuous phase transitions [An¬ 
dreev 1982] and performs an analytic expansion of /p(p) 
in p ]Cabrera 1964, Cabrera and Garcia 1982] (and trun¬ 
cate after a quadratic term f 2 P^), then a similar pro¬ 
cedure leads d = 2, which is often referred to as the 
“mean-field” value. This same value can be produced by 
quenched impurities, as shown explicitly for the equiv¬ 
alent commensurate-incommensurate transition by Kar- 
dar and Nelson [1985]. 

There are some other noteworthy results for the 
smooth edge. As the facet roughening temperature is 
approached from below, the facet radius shrinks like 
exp[—7r^TR/4{21n2(rR — T)}^/^] [Jayaprakash et al. 
1983], in striking contrast to predictions by mean-field 
theory. The previous discussion implicitly assumes that 
the path along x for which d = 3/2 in Eq. (33) is nor¬ 


mal to the facet edge. By mapping the crystal surface 
onto the asymmetric 6-vertex model, using its exact so¬ 
lution [Yang 1967, Sutherland et al. 1967], and employ¬ 
ing the Bethe Ansatz to expand the free energy close to 
the facet edge, Dahmen et al. [1998] find that d = 3/2 
holds for any direction of approach along the rounded 
surface toward the edge, except along the tangential di¬ 
rection (the contour that is tangent to the facet edge at 
the point of contact Xq- In that special direction, they 
find the new critical exponent dy = 3 (where the sub¬ 
script y indicates the direction perpendicular to the edge 
normal, x. Akutsu et al. [1998] (also Akutsu and Akutsu 
[2006]) confirmed that this exact result was universally 
true for the Gruber-Mullins-Prokrovsky-Talapov free- 
energy expansion. (The Pokrovsky-Talapov argument 
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was for the equivalent commensurate-incommensurate 
transition.) They also present numerical confirmation 
using their transfer-matrix method based othe product- 
wave-function renormalization group (PWFRG) [Nishino 
and Okunishi 1995; Okunishi et al. 1999]. Observing -dy 
experimentally will clearly be difficult, perhaps impossi¬ 
ble; the nature and breadth of crossover to this unique 
behavior has not, to the best of my knowledge, been pub¬ 
lished. A third result is that there is a jump (for T < Tr) 
in the curvature of the rounded part near the facet edge 
that has a universal value [Akutsu et al. 1988; Sato and 
Akutsu 1995], distinct from the universal curvature jump 
of the ECS at Tr [Jayaprakash et al. 1983]] 


6.2. Experiments on Lead 


Noteworthy initial experimental tests of d = 3/2 
include direct measurements of the shape of equili¬ 
brated crystals of "^He [Carmi et al. 1987] and Pb 
[Rottman...Metois 1984]. As in most measurements of 
critical phenomena, but even harder here, is the iden¬ 
tification of the critical point, in this case the value of 
xq at which rounding begins. Furthermore, as is evident 
from Eq. (33), there are corrections to scaling, so that 


the “pure” exponent 3/2 is seen only near the edge and 
a larger effective exponent will be found farther from the 
edge. For crystals as large as a few mm at tempera¬ 
tures in the range 0.7K-1.1K, "^He d = 1.55 ± 0.06 was 
found, agreeing excellently with the Pokrovsky-Talapov 
exponent. The early measurements near the close-packed 
(111) facets of Pb crystallites, at least two orders of mag¬ 
nitude smaller, were at least consistent with 3/2, stated 
conservatively as d = 1.60 ± 0.15 after extensive analy¬ 
sis. Saenz and Garcia [1985] proposed that in Eq. (28) 
there can be a quadratic term, say f 2 P^ (but neglect the 
possibility of a quartic term). Carrying out the Legen¬ 
dre transformation then yields an expression with both 
X — B and {x — B + terms which they claim 

will lead to effective values of d between 3/2 and 2. This 
approach provided a competing model for experimental¬ 
ists to consider but in the end seems to have produced 
little fruit. 

As seen in Fig. STM allows detailed measurement 
of micron-size crystal height contours and profiles at fixed 
azimuthal angles. By using STM to locate the initial 
step down from the facet, first done by Surnev et al. 
[1998] for supported Pb crystallites, xq can be located 
independently and precisely. However, from the 1984 
Heyraud-Metois experiment [Rottman...Metois 1984] it 
took almost two decades until the Bonzel group could be 
fully confirm the d = 3/2 behavior for the smooth edges 
of Pb(lll) in a painstaking study [Nowicki et al. 2002b]. 
There were a number of noteworthy challenges. While 
the close-packed 2D network of spheres has 6-fold sym¬ 
metry, the top layer of a (111) facet of an fee crystal (or 
of an (0001) facet of an hep crystal) has only 3-fold sym¬ 
metry due to the symmetry-breaking role of the second- 



FIG. 18: (a) Micron-size lead crystal (supported on Ru) im¬ 
aged with a variable-temperature STM at T = 95° C. Anneal¬ 
ing at T = 95° C for 20 hours allowed it to obtain its stable, 
regular shape. Lines marked A and B indicate location of 
proHles. Profile A crosses a (0 0 l)-side facet, while profile B 
a (1 1 l)-side facet, (b) 770 nmx770 nm section of the top 
part of a Pb-crystal. The insert shows a 5.3 nmx5.3 nm area 
of the top facet, confirming its (1 1 l)-orientation. Both the 
main image and the insert were obtained at T = 110°C. From 
Thiirmer et al. [2003]. 


layer. There are two dense straight step edges, called A 
and B, with {100} and {111} microfacets, respectively. 
In contrast to noble metals, for Pb there is a sizeable (of 
order 10%) difference between their energies. Even more 
significant, when a large range of polar angles is used in 
the fitting, is the presence of small (compared to (111)) 
{112} facets for equilibration below 325K. Due to the 
high atomic mobility of Pb that can lead to the forma¬ 
tion of surface irregularities, Surnev et al. [1998] worked 
close to room temperature. One then finds strong (3-fold) 
variation of d with azimuthal angle, with § oscillating 
between 1.4 and 1.7. With a higher annealing tempera¬ 
ture of 383K Nowicki et al. [2002b] report the azimuthal 
averaged value § = 1.487 (but still with sizeable oscilla¬ 
tions of about ±0.1); in a slightly early short report, they 
[Nowicki et al. 2002a] give a value d = 1.47 for annealing 
at room temperature. Their attention shifted to deduc¬ 
ing the strength of step-step repulsions by measuring g 
[Nowicki et al. 2003, Bonzel and Nowicki 2004]. In the 
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most recent review of the ECS of Pb, Bonzel et al. [2007] 
rather tersely report that the Pokrovsky-Talapov value 
of 3/2 for characterizes the shape near the (111) facet 
and that imaging at elevated temperature is essential to 
get this result; most of their article relates to comparison 
of measured and theoretically calculated strengths of the 
step-step interactions. 

Few other systems have been investigated in such de¬ 
tail. Using scanning electron microscopy (SEM) Metois 
and Heyraud [1987] considered In, which has a tetrago¬ 
nal structure, near a (111) facet. They analyzed the re¬ 
sulting photographs from two different crystals, viewed 
along two directions. For polar angles 0° < d < ^ 5° 
they find « 2 while for 5° < 0 < 15° determine 
0 « 1.61, concluding that in this window 0 = 1.60±0.10; 
the two ranges have notably different values of xq. This 
group [Bermond et al, 1998] also studied Si, equilibrated 
at 900°C, near a (111) facet. Many profiles were mea¬ 
sured along a high-symmetry (111) zone of samples with 
various diameters of order a few /rm, over the range 
3° < 0 < 17°. The results are consistent with d = 3/2, 
with an uncertainty estimated at 6%. Finally Gladic et 
al. [2002] studied large (several mm.) spherical cuprous 
selenide (Cu 2 -a;Se) single crystals near a (111) facet. 
Study in this context of metal chalcogenide superionic 
conductors began some dozen years ago because, other 
than "^He, they are the only materials having sub-cm. 
size crystals that have an ECS form that can be grown 
on a practical time scale (viz. over several days) because 
their high ionic and electronic conductivity enable fast 
bulk atomic transport. For 14.0° < 9 < 17.1° Gladic 
et al. [2002] find = 1.499 ± 0.003. (They also report 
that farther from the facet d « 2.5, consistent with the 
Andreev mean field scenario.) 


6.3. Summary of Highlights of Novel Approach to 
Behavior Near Smooth Edges 


Digressing somewhat, we note that Ferrari, Prahofer, 
and Spohn [2004] (FPS) found novel static scaling be¬ 
havior of the equilibrium fluctuations of an atomic ledge 
bordering a crystalline facet surrounded by rough regions 
of the ECS in their examination of a 3D Ising corner 
(Fig. 19). This boundary edge might be viewed as a 
“shoreline” since it is the edge of an island-like region-the 
crystal facet-surrounded by a “sea” of steps [Pimpinelli 
et al. 2005]. 

FPS assume that there are no interactions between 
steps other than entropic, and accordingly map the step 
configurations can be mapped to the world lines of free 
spinless fermions, as in treatments of vicinal surfaces 
[Jayaprakash et al. 1984]. However, there is the key new 
feature that the step number operator is weighted by the 
step number, along with a Lagrange multiplier asso¬ 
ciated with volume conservation of the crystallite. The 
asymmetry of this term is leads to the novel behavior 
they find. They then derive an exact result for the step 



FIG. 19: Simple-cubic crystal corner viewed from the {111} 
direction. From Ferrari et al. [2004]. 


density and find that, near the shoreline, 

lim = —x(Ai(a;))^ -I- (Ai'(a;))^, (34) 

A—>-oo 

where px is the step density (for the particular value of 

A). 

The presence of the Airy function Ai results from the 
asymmetric potential implicit in Up and preordains ex¬ 
ponents involving 1/3. The variance of the wandering of 
the shoreline, the top fermionic world line in Fig. and 
denoted by b, is given by 

Var[6A(t) - bx{0)] = X^/^g{X~'^/^t) (35) 

where t is the fermionic “time” along the step; g{s) ~ 2jsj 
for small s (diffusive meandering) and ^ 1.6264 — 2/s^ 



FIG. 20: Magnified detail of the steps near the facet edge, 
from Ferrari and Spohn [2003]; b) Snapshot of computed con¬ 
figurations of the top steps (those near a facet at the flattened 
side portion of a cylinder) for a terrace-step-kink (TSK) model 
with volume constraint. From Ferrari et al. [2004]. 
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FIG. 21: Summary of experimental results for vicinal Si (111) 
surface: Q denotes the temperature at which faceting begins 
for surfaces misoriented towards the ( 110 ] direction, x the 
faceting temperatures for surfaces misoriented towards the [ 
112 ], and □ the temperatures at which the step structure of 
surfaces misoriented towards the [112] direction ch ang e. The 
dashed line displays a fit of the 110] data to Eq. (40l. The 
dotted lines show how a 4° sample phase-separates into the 
states denoted by • as it is further cooled. From Bartelt et 
al. [1989]. 


for large s. They then 1.202... is Apery’s constant and 
N is the number of atoms in the crystal. They find 

Yar[bi{eT + x)-bi{eT)] ^ , (36) 

where A = {l/2)b'^. This leads to their central result 
that the width w ^ in contrast to the scaling 
of an isolated step or the boundary of a single-layer island 
and to the ln£ scaling of a step on a vicinal surface, i.e. 
in a step train. Furthermore, the fluctuations are non- 
Gaussian. They also show that near the shoreline, the 
deviation of the equilibrium crystal shape from the facet 
plane takes on the Pokrovsky-Talapov [1979, 1984[ form 
with d = 3/2. 

From this seminal work we could derive the dynamic 
exponents associated with this novel scaling and mea¬ 
sure them with STM [Pimpinelli et al. 2005, Degawa et 
al. 2006, 2007], as reviewed in Einstein and Pimpinelli 
[2014]. 


7. SHARP EDGES AND FIRST-ORDER 
TRANSITIONS—EXAMPLES AND ISSUES 

7.1. Sharp Edges Induced by Facet Reconstruction 

Si near the (111) plane offers an easily understood 
entree into sharp edges [Phaneuf and Williams 1987, 
Bartelt et al. 1989]. As Si is cooled from high tempera¬ 
tures, the ( 111 ) plane in the “(1 x 1 )” phase reconstructs 


Phase Diagram Wulff Plo* 



FIG. 22: Wulff plots illustrating the effect of a reconstruc¬ 
tive transition on the EGS, and corresponding temperature- 
[mis]orientation phase diagrams. The solid curves represent 
the EGS with an unreconstructed [“( 1 x 1)”] facet, while the 
dashed curves give the EGS with a reconstructed facet. As 
temperature decreases, the free energy of the reconstructed 
facet, relative to that of the unreconstructed facet, decreases. 
Below the transition temperature Tc (called T 7 in the text), 
the two shapes intersect, giving a “net” EGS that is the in¬ 
ner envelope of the two. The phase diagram shows regions 
where all orientations tan 0 (or m) are allowed for the unre¬ 
constructed crystal [“(1 x 1 )”], regions of phase separation 
(labeled “coex.”), and regions where the reconstruction (la¬ 
beled “rec.”) is allowed for ranges of orientation. The relative 
size of the reconstructed and unreconstructed facets depends 
on the free energy to create a step on the reconstructed ( 111 ) 
face, compared to its unreconstructed counterpart: (a) the 
behavior for extremely large energy to create steps on the 
(7 X 7) terrace and (b) a smaller such energy. Solid circles 
mark the sharp edge at the temperature at which the crystal 
shapes cross. Crosses show the intersection of the facet and 
the curved part (i.e., the smooth edge) of the crystal shape 
for the reconstructed phase. From Bartelt et al. [1989]. 


into a (7 X 7) pattern [Schlier and Farnsworth 1959] 
around 850°C, to be denoted T 7 to distinguish it clearly 
from the other subscripted temperatures. (The notation 
“(1 X 1)” is intended to convey the idea that this phase 
differs considerably from a perfect ( 111 ) cleavage plane 
but has no super lattice periodicity.) For comparison, the 
melting temperature of Si is ^ 1420° C, and the Tr is 


estimated to be somewhat higher. As shown in Fig. 21 


above T 7 surfaces of all orientations are allowed and are 
unreconstructed. At T/ a surface in the (111) direction 
reconstructs but all other orientations are allow and are 
unreconstructed. Below Tj, for surfaces misoriented to¬ 
ward [ 112 ] remain stable during cooling (although the 
step structure changes). On the other hand, on surfaces 
misoriented toward [IlO] and [112] the temperature at 
which the (7 x 7) occurs decreases with increasing mis- 
orientation angle rh. Furthermore, just as the (7 x 7) 
appears, the surface begins to separate into two phases, 
one a perfectly oriented (7 x 7) plane, m = 0 and the 
second an unreconstructed phase with a misorientation 
greater than that at higher temperature. As tempera¬ 
ture further decreases, the misorientation of the unrecon¬ 
structed phase increases. Figure depicts this scenario 
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with solid circles and dotted lines for a 4° misoriented 
sample at 840° C. This behavior translates into the for¬ 
mation of a sharp edge on the ECS between a flat (7 x 
7) line and a rounded “(1 x 1)” curve. 

To explain this behavior, one co-plots the ECS for the 
two phases, as in Fig. 22 [Bartelt et al. 1989]. The free 
energy to create a step is greater in the (7 x 7) than in the 
“(1 X 1)” phase. In the top panels (a), the step energy 
for the (7 x 7) is taken as infinite (i.e. much larger than 
that of the “(1 x 1)” phase) so its ECS never rounds. 
At T^ (Tc on the figure) the free energies per area /o 
of the two facets are equal, call them /y and /i, with 
associated energies uy and ui and entropies sy and si for 
the (7 X 7) and “(1 x 1)” phases, respectively, near Ty. 
Then Ty = (uy — uy)/(si — Sy) and, assuming the internal 
energies and entropies are insensitive to temperature. 


/i -/7 = (Ty-T)(si-sy). (37) 

Since si > sy because the (7 x 7) phase is so highly 
ordered, we find that /i — /y > 0 below Ty, as illustrated 
in Fig. Making connection to thermodynamics, we 
identify 


L 

TV 




(38) 


where L is the latent heat of the first-order reconstruction 
transition. 

Corresponding to the minimum of a free energy as dis¬ 
cussed earlier, the ECS of the system will be the inner 
envelope of the dashed and solid traces: a flat (7 x 7) 
facet along the dashed line out to the point of intersec¬ 
tion, the sharp edge, beyond which it is “(1 x 1)” with 
continuously varying orientation. If one tries to construct 
a surface with a smaller misorientation, it will phase sep¬ 
arate into flat (7 x 7) regions and vicinal unreconstructed 
regions with the orientation at the curved (rough) side of 


the sharp edge. Cf. Fig. 23 


Using the leading term in Eq. (32) or (33) we can es¬ 


timate the slope of the coexisting vicinal region and its 
dependence on temperature [3]: First we locate the sharp 
edge (recognizing /o as /i and zq as zi ) by noting 

Z7 = zi - 2(A/g)^/^(a; - xo)^/^ (39) 

(Ty-T)As « (/i-/y)r = A3/V'/"(^-a;o)"/" 

Since the slope m there is the 

temperature dependence of the slope is 


=‘'1) 


1/3 


1--) 

Tr) 


1/3 


(40) 


If the step free energy of the reconstructed phase were 
only modestly greater than that of the “(1 x I)”, then, as 
shown in the second panel in Fig. the previous high- 
T behavior obtains only down to the temperature Ty at 
which the “(1 x I)” curve intersects the (7 x 7) curve at 
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FIG. 23: Microscopic view of what happens to a misoriented 
surface in Fig. as temperature decreases, (a) At high tem¬ 
perature, the Si(Ul) vicinal surface is a single, uniform phase. 
Initial terrace widths t are typically a few nm., as determined 
by the net angle of miscut oo (i.e., do), and the step-height 
h, which is one interplanar spacing (~ 0.31 nm). (b) Below 
the (7x7) reconstruction temperature (~ 850°C) the steps 
cluster to form a new surface of misorientation angle a(T) 
(i.e., 6^). A facet of (111) orientation with (7x7) reconstruc¬ 
tion forms simultaneously. The width of the (111) facet, £, 
is larger than the experimentally resolvable width of 500 A. 
(c) Well below the transition, the step separation reaches a 
minimum distance, tmin ~ 1 nm. No further narrowing oc¬ 
curs, perhaps because surface diffusion is too slow T < 600°C. 
From Williams and Bartelt [1989]. 


its [smooth] edge. For T < Ti the sharp edge associated 
with the interior of the curves is between a misoriented 
“(1 X 1)” phase and a differently misoriented (7 x 7) 
phase, so that it is these two which coexist. All orienta¬ 
tions with smaller misorientation angles than this (7 x 
7) plane are also allowed, so that the forbidden or coex¬ 
istence regime has the depicted slivered crescent shape. 
Some other, but physically improbable, scenarios are also 
discussed by Bartelt et al. [1989]. Phaneuf and Williams 
[1987] show (their Fig. 3) the LEED-beam splitting for 
a surface misoriented by 6.4° is oc (Ty — T)^/^ once the 
surface is cooled below the temperature (which is < Ty) 
when this orientation becomes unstable to phase sepa¬ 
ration; however, by changing the range of fitting, they 
could also obtain agreement with (Ty — T)^/^, i.e. = 
2. With high-resolution FEED Hwang et al. [1989] con- 
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elude that the exponent /3 = {'& — = 0.33 ± 0.05 

(i.e., that d = 3/2. The result does depend somewhat 
on what thermal range is used in the fit, but they can 
decisively rule out the mean-field value ■§ = 2. Williams 
et al. [1993] give a more general discussion of vicinal Si, 
with treatment of azimuthal in addition to polar mis- 
orientations. In contrast, synchrotron x-ray scattering 
experiments by Noh et al. [1991, 1993] report the much 
larger d = 2.3 ipj. However, subsequent synchrotron 
x-ray scattering experiments by Held et al. [1995] obtain 
a decent fit of data with '0 = 3/2 and a best fit with 

= 1.75 (i.e., P = 0.43 ± 0.07). (They also report that 
above 1159K, the surface exists as a single, logarithmi¬ 
cally rough phase.) The origin of the curious value of i? 
in the Noh et al. experiments is not clear. It would be 
possible to attribute the behavior to impurities, but there 
is no evidence to support this excuse, and indeed for the 
analogous behavior near the reconstructing (331) facet 
of Si (but perhaps a different sample), Noh et al. [1995] 
found 'd = 1.47 ±0.1. It is worth noting that extracting 
information from x-ray scattering from vicinal surfaces 
requires great sophistication (cf. the extensive discussion 
in Dietrich and Haase [1995]) and attention to the size of 
the coherence length relative to the size of the scattering 
region [Ocko, 2014], as for other diffraction experiments. 

Similar effects to reconstruction (viz. the change in /o) 
could be caused by adsorption of impurities on the facet 
[Cahn 1982]. Some examples are given in a review by So- 
morjai and Van Hove [1989]. In small crystals of dilute 
Pb-Bi-Ni alloys, co-segregation of Bi-Ni to the surface 
has a similar effect of reversibly) changing the crystal 
structure to form {112} and {110} facets [Cheng and 
Wynblatt, 1994]. There is no attempt to scrutinize the 
ECS to extract an estimate of d. Meltzman et al. [2011] 
considered the ECS of Ni on a sapphire support, noting 
that, unlike most fee crystals, it exhibits a faceted shape 
even with few or no impurities, viz. with {111}, {100}, 
and {110} facets; {135} and {138} emerged at low oxy¬ 
gen pressure and additionally {012} and {013} at higher 
pressure. 


The phase diagram of Pt(OOl), shown in Fig. 24 and 
studied by Yoon et al. [1994] using synchrotron x-ray 
scattering, at first seems similar to that of Si near (111) 
[Song et al. 1994, 1995], albeit with more intricate magic 
phases with azimuthal rotations at lower temperatures, 
stabilized by near commensurability of the period of their 
reconstruction and the separation of their constituent 
steps. In the temperature-misorientation (surface slope) 
phase diagram, shown in Fig. 24 the (001) facet under¬ 
goes a hexagonal reconstruction at Tg = 1820 K (well be¬ 
low the bulk melting temperature of 2045K). For samples 
misoriented from the (100) direction (which are stable at 
high temperature), there is coexistence between flat re¬ 
constructed Pt(OOl) and a rough phase more highly mis¬ 
oriented than it was at high temperature, with a misori- 
entation that increases as temperature decreases. How¬ 
ever, they hnd /3 = 0.49 ± 0.05, or r? = 1.96, consistent 
with mean field and inconsistent with /3 = l/3orr? = 3/2 



FIG. 24: Orientational phase diagram of vicinal Pt (001) mis¬ 
oriented toward the [110] direction. Single-phase regions are 
hatched, and two-phase coexistence regions are unhatched. 
Solid lines are boundaries between two phases. Dashed lines 
mark triple points. Open circles show misorientation angles 
measured for a sample miscut by 1.4° towards the [110] direc¬ 
tion, while solid circles show tilt angles measured for a sample 
miscut by 3.0°. From Yoon et al. [1994]. 


of Pokrovsky-Talapov. The source of this mean-held ex¬ 
ponent is that in this case the (001) orientation is rough 
above Tg. Hence, in Eq. {28) B vanishes, leaving the ex¬ 


pansion appropriate to rough orientations . Proceeding 
as before, Eq. (32 becomes 


fp{p) = fo + fiv) = fo- 


(41) 


where the result for f{r]) is reached by proceeding as 
before to reach the modiheation of Eq. ( [32| ). Thus, there 
is no smooth edge take-off point (no shoreline) in the 
equivalent of Fig. [^ and one hnds the reported exponent 

near 2. 

The effect of reactive and nonreactive gases metal cat¬ 
alysts has long been of interest [Flytzani-Stephanopoulos 
and Schmidt, 1979]. Various groups investigated 
adsorbate-induced faceting. Walko and Robinson [2001] 
considered the oxygen-induced faceting of Cu(115) into 
O/Cu(104) facets, using Wulff constructions to ex¬ 
plain their observations. They found three temperature 
regimes with qualitatively different faceting processes. 
Szczepkowicz et al. [2005] studied the formation of {211} 
facets by depositing oxygen and paladium on tungsten, 
both on (111) facets and on soherical crystals. While the 
shape of the facets is different for flat and curved surfaces, 
the distance between parallel facet edges is comparable, 
although the area of a typical facet on a curved surface 
is an order of magnitude greater. There is considerable 
information about facet sizes, width of the facet-size dis¬ 
tribution, and surface rms roughness. 

For 2D structures on surfaces, edge decoration can 
change the shape of the islands. A well-documented ex¬ 
ample is Pt on Pt(lll). As little as 10“^ ML of CO pro¬ 
duces a 60° rotation of the triangular islands by change 
the balance of the edge free energies of the two differ¬ 
ent kinds of steps forming the island periphery [Kalff et 
al. 1998]. Stasevich et al. [2009] showed how decoration 










20 


of single-layer Ag islands on Ag(lll) by a single-strand 
“necklace” of Ceo dramatically changes the shape from 
hexagonal to circular. With lattice-gas modeling com¬ 
bined with STM measurements, they could estimate the 
strength of Ceo-Ag and Ceo-Ceo attractions. Generaliza¬ 
tions to decoration on systems with other symmetries is 
also discussed. 


8. GOLD-PROTOTYPE OR ANOMALY OF 
ATTRACTIVE STEP-STEP INTERACTION? 


Much as "^He and Pb are the prototypical materials 
with smooth edges, Au is perhaps the prime example of 
a surface with sharp edges, around (111) and (100) facets. 
(Cf., e.g., Wortis [1988]). Care must be taken to insure 
that the surface is not contaminated by atoms (typically 
C) from the supporting substrate [Wang and Wynblatt 
1988]. (See similar comments by Handwerker et al. [1988] 
for ceramics, which have a rich set of ECS possibilities.) 
To describe these systems phenomenologically, the pro¬ 
jected free energy expansion in Eq. 0 requires a negative 
term to generate a region with negative curvature, as in 
Fig.0 so that the two orientations joined by the Maxwell 
double-tangent construction correspond to the two sides 
of the sharp edge. Thus, for sharp edges around facets, 
the more-left minimum must be in the high-symmetry 
facet direction. 

In a mean field based approach, Wang and Wynblatt 
[1988] included a negative quadratic term, with ques¬ 
tionable physical basis. Emundts et al. [2001] instead 
took the step-step interaction to be attractive {g < 0) in 
Eq. §. Then proceeding as above they find 





(42) 


where Pc is the tangent of the the facet contact angle. 
Note that both the shift in the facet edge from B and the 
contact slope increase with \g\/c. Emundts et al. [2001] 
obtain estimates of the key energy parameters in the ex¬ 
pansion for the sharp edges of both the (111) and (100) 
facets. They also investigate whether it is the lowering 
of the facet free energy /o that brings about the sharp 
edges, in the manner of the case of Si(lll) discussed 
above. After reporting the presence of standard step-step 
repulsions (leading to narrowing of the TWD) in experi¬ 
ments on flame-annealed gold, Shimoni et al. [2000] then 
attribute to some effective long-range attraction—with 
undetermined dependence on ^—the (non-equilibrium) 
movement of single steps toward step bunches whose 
steps are oriented along the high-symmetry (110). 

Is it possible to find a generic long-range attractive 
A£~'^ step-step interaction {A < 0) for metals and ele¬ 
mental semiconductors (where there is no electrostatic 
attraction between oppositely charged atoms)? Several 
theoretical attempts have only been able to find such at¬ 
tractions when there is significant alternation between 
“even” and “odd” layered steps. Redfield and Zangwill 


consider whether surface relaxation can produce such an 
attraction, pointing out a flaw in an earlier analysis as¬ 
suming a rigid relaxation by noting that for large step 
separations, the relaxation must return to its value for 
the terrace orientation. Since atomic displacements fall 
off inversely with distance from a step, the contribution 
to the step interaction can at most go like and tend 
to mitigate the combined entropic and elastic repulsion. 
They argue that this nonlinear effect is likely to be small, 
at least for metals. It is conceivable that on an elastically 
highly anisotropic surface, the elastic interaction might 
not be repulsive in special directions, though I am not 
aware of any concrete examples. 

By observing that the elastic field mediating the in¬ 
teraction between steps is that of a dipole applied on 
a stepped rather than on a flat surface, Kukta et al. 
[2002] deduce a correction to the behavior of the 
Marchenko-Parshin [1980] formula that scales as ln£. 
Using what was then a state-of-the-art semiempirical po¬ 
tential, the embedded atom method (EAM)[Daw et al. 
1993], they find that this can lead to attractive interac¬ 
tions at intermediate values of £. However, their “rough¬ 
ness correction” term exists only when the two steps have 
unlike orientations (i.e., one up and one down, such as on 
opposite ends of a monolayer island or pit). For the like- 
oriented steps of a vicinal surface or near a facet edge, 
the correction term vanishes. The oft-cited paper then 
invokes 3-step interactions, which are said to have the 
same size as the correction term, as a way to achieve at¬ 
tractive interactions. Although the authors discuss how 
this idea relates to the interaction between an isolated 
step and a step bunch, they do not provide the explicit 
form of the threefold interaction; their promise that it 
will be “presented elsewhere” has not, to the best of my 
searching, ever been fulfilled. Prevot and Croset [2004] 
revisited elastic interactions between steps on victuals 
and found that with a buried-dipole model (rather than 
the surface-dipole picture of Marchenko and Parshin), 
they could achieve “remarkable agreement” with molec¬ 
ular dynamics simulations for victuals to Cu and Pt (001) 
and (111), for which data is fit by 
The tabulated values of E^^ indeed agree well with their 
computed results for their improved elastic model, which 
includes the strong dependence of the interaction energy 
on the force direction. While there is barely any discus¬ 
sion of i? 3 , plots of the interaction are always repulsive. 
Hecquet [2008] finds that surface stress modifies the step- 
step interaction compared to the Marchenko-Parshin re¬ 
sult, enhancing the prefactor of nearly threefold for 
Au(OOl); again, there is no mention of attractive inter¬ 
actions over any range of step separations. 

In pursuit of a strictly attractive step interaction 
to explain the results of Shimoni et al. [2000], Wang et 
al. [2007] developed a model based on the SSH model 
[Su et al. 1979] of polyacetylene (the original model ex¬ 
tended to include electron-electron interaction), focusing 
on the dimerized atom rows of the (2 x 1) reconstruc¬ 
tion of Si(OOl). The model produces an attractive cor- 
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rection term to the formula derived by Alerhand et al. 
[1988] for interactions between steps on Si(OOl), where 
there is ABAB alternation of (2 x 1 and 1x2) recon¬ 
structions on neighboring terraces joined at single-height 
steps. For this type of surface, the correction has lit¬ 
tle significance, being dwarfed by the logarithmic repul¬ 
sion. It also does not occur for vicinals to high-symmetry 
facets of metals. However, for surfaces such as Au(llO) 
with its missing row morphology [Copel and Gustafsson, 
1986] or adsorbed systems with atomic rows, the row can 
undergo a Peierls [1955] distortion that leads to an anal¬ 
ogous dimerization and an l~'^ attraction. There have 
been no tests of these unsung predictions by electronic 
structure computation. 

Returning to gold, applications of the glue potential (a 
semiempirical potential rather similar to EAM), Ercolessi 
et al. [1987] were able to account for reconstructions of 
various gold facets, supporting that the sharp edges on 
the ECS are due to the model used for Si(lll) rather than 
attractive step interactions. Studies by this group found 
no real evidence for attractive step interactions [Tosatti, 
2014]. 

In an authoritative review a decade ago, Bonzel 
[2003]—the expert in the field who has devoted the 
most sustained interest in ECS experiments on elemen¬ 
tal systems—concluded that it was not possible to de¬ 
cide whether the surface reconstruction model or attrac¬ 
tive interactions was more likely to prove correct. In my 
view, mindful of Ockham’s razor, the former seems far 
more plausible, particularly if the assumed attractive in¬ 
teraction has the form. 

The phase diagram of surfaces vicinal to Si(113) 
presents an intriguing variant of that vicinal to Si(lll). 
There is again a coexistence regime between the (113) ori¬ 
entation and progressively more highly misoriented vici¬ 
nals as temperature is reduced below a threshold temper¬ 
ature T(, associated with a first-order transition. How¬ 
ever, for higher temperatures T > Tj there is a continu¬ 
ous transition, in contrast to the behavior on (111) sur¬ 
faces for T > Tt. Thus, Song and Mochrie [1995] identify 
the point along (113) at which coexistence vanishes, i.e. 
Tt, as a tricritical point, the first such point seen in a 
misorientation phase diagram. To explain this behavior 
Song and Mochrie invoke a mean-field Landau-theory ar¬ 
gument in which the cubic term in p is proportional to 
{T — Tt), so negative for T < Tt, with a positive quar- 
tic term. Of course, this produces the observed generic 
behavior, but the exponent /3 is measured as 0.42± 0.10 
rather than the mean-field value 1. Furthermore, the 
shape of the phase diagram differs from the mean field 
prediction and the amplitude of the surface stiffness be¬ 
low Tt is larger than above it, the opposite of what hap¬ 
pens in mean held. Thus, it is not clear in detail what 
the interactions actually are, let alone how an attractive 
interaction might arise physically. 


9. WELL-ESTABLISHED ATTRACTIVE 
STEP-STEP INTERACTIONS OTHER THAN 

For neutral crystals, there are two ways to easily ob¬ 
tain interactions that are attractive for some values of 
In neither case are the interactions monotonic long-range. 
The hrst is short-range local effects due to chemical prop¬ 
erties of proximate steps while the other is the indirect 
Friedel-like interaction. 


9.1. Atomic-Range Attractions 


At very small step separations the long-range the long- 
range monotonic behavior is expected to break down 
and depend strongly on the local geometry and chem¬ 
istry. Interactions between atoms near step edges are 
typically direct and so stronger than interactions medi¬ 
ated by substrate elastic helds or indirect electronic ef¬ 
fects (see below). We saw earlier that a higher-order 
term arises at intermediate separations [Najafabadi and 
Srolovitz 1994], and further such terms should also ap¬ 
pear with decreasing On TaC(910) [vicinal to (001) 
and miscut toward the [010] direction], Zuo et al. [2001] 
explained step bunching using a weak £-3[±0'5] attrac¬ 
tion in addition to the repulsion. [The double-height 
steps are electrically neutral.] Density-functional theory 
(DFT) studies were subsequently performed for this sys¬ 
tem by Shenoy and Ciobanu [2003]. Similarly, Yamamoto 
et al. [2010] used an attractive dipole-quadrupole in¬ 
teraction to explain anomalous decay of multilayer holes 
on SrTiO3(001). 

More interesting than such generic effects are attrac¬ 
tions that occur at very short step separations for special 
situations. A good example is Ciobanu et al. [2003], who 
find an attraction at the shortest separation due to the 
cancellation of force monopoles of two adjacent steps on 
vicinal Si(113) at that value of 

As alluded to above, most of our understanding of the 
role of step interactions comes from the mapping of 
classical step configurations in 2D to the world lines of 
spinless fermions in ID. Unlike fermions, however, steps 
can touch (thereby forming double-height steps), just not 
cross. Such behavior is even more likely for vicinal fee 
or bcc (001) surfaces, where the shortest possible “ter¬ 
race”, some fraction of a lateral nearest-neighbor spac¬ 
ing, amounts to touching fermions when successive lay¬ 
ers of the crystal are described with simple-cubic rather 
than layer-by-layer laterally offset coordinates. Sathiya- 
narayanan et al. [2009] investigated some systematics of 
step touching, adopting a model in which touching steps 
on a vicinal cost an energy ct- Note that et = oo recoups 
the standard fermion model. For simplicity, the short 
study concentrates on the “free fermion” case A = 0, i.e. 
g = 2. (Cf. Eq. (25).) Even for et = 0 there is an effec¬ 
tive attraction, i.e. g <2, since the possibility of touching 
broadens the TWD. This broadening is even more pro¬ 
nounced for et < 0. In other words, such short-range 
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effects can appear, for a particular system, to contribute 
a long-range attraction. Closer examination shows that 
this attraction is a finite-size effect that fades away for 
large values of (£). In our limited study, we found that 
fits of simulated data to the GWD expression could be 
well described by the following finite-size scaling form, 
with the indicated three fitting parameters: 


Peff = 2-(0.9±0.1)(^)-°'29±oo7exp[-(3.3±0.2)et/fcBT]. 

^ (43) 

While Eq. (431 suggests that making the step touch¬ 
ing more attractive (decreasing ct) could decrease /3eff 
without limit, instabilities begin to develop, as expected 
since Lassig [1996] showed that for A < —1/4, i.e. 
A < —{kBT)‘^/AI3, a vicinal surface becomes unstable 
(to collapse to step bunches). Correspondingly, the low¬ 
est value tabulated in Sathiyanarayanan et al. [2009] is 
etlkBT =-0.2. 

To distinguish true long-range {£~^) attractions on vic¬ 
inal surfaces requires measurements of several different 
vicinalities (i.e. values of {£)). Likewise, in analyses of 
ECS data, consideration of crystallites of different sizes 
would seem necessary. Wortis [1988] had noted the im¬ 
portance of size dependence in other contexts. 

Along this theme, an instructive specihc case is the 
“sticky-step” or, more formally, the p-RSOS [restricted 
solid-on-solid with point-contact attractions between 
steps] model explored in detail by Akutsu [2011] using the 
product wavefunction renormalization group (PWFRG) 
method, calculating essentially the ECS (see Fig. 
and related properties. Steps are zig-zag rather than 
straight as in the preceding Sathiyanarayanan model, so 
her “stickiness” parameter eint is similar but not identi¬ 
cal to et. She finds that in some temperature regimes, 
non-universal non-Pokrovsky Talapov values of £} occur. 
Specifically, let T/^iii(ei„t/e, (/o) and T/^ooi(emt/e, (/>o) be 
the highest temperature at which a first-order phase tran¬ 
sition (sharp edge) occurs for the (111) and (001) facets, 
respectively, where indicates the position along the 

ECS. Note Tfiii{et/e,(j)o) = (0.3610 ± 0.0005)e/fcB > 
T/,iii(e™t/e, </o) = (0.3585 ± 0.0007)e/fcB. For fcsT/e = 
0.37, so T > T/_iii(—0.5,7r/4), Akutsu recovers 
Pokrovsky-Talapov values for £) and £}t, but for kBT/e = 
0.36 (shown in Fig. 26), so T/ 4 ii(—0.5,7r/4) > T > 





FIG. 26: Profiles in the diagonal direction of the surface 
in Fig. |25| still at ksTjei = 0.3. Broken lines represent 
metastable lines, (a) ksT/ei = 0.36, tint jet = —0.5, on a 
very fine length scale. The edge of the (111) facet is denoted 
by Xq. (b) The original RSOS model {eint = 0) on a much 
coarser scale. On this scale (and on an intermediate scale not 
included here), the profile profiles are flat until the edge. On 
the intermediate scale, the region beyond Xg is starts deviat¬ 
ing rather smoothly for ksT/ei = 0.35 but looks straight for 
ksT/ei = 0.36 and 0.37. See text and source. From Akutsu 
[ 2011 ]. 


Tf, ooi(—0.5,7r/4), the values are very different: £) = 
1.98±0.03 and £)t = 3.96±0.08, more like mean field. For 
4>o = 0 (tilting toward the (100) direction) only standard 
Pokrovsky-Talapov exponents are found. Upon closer ex¬ 
amination with Monte Carlo simulations, Akutsu finds 
large step bunches for T < T/qoo but step droplets for 
Tf ,001 < T < Tf,iii. The details are beyond the scope 
of this review, but eventually Akutsu deduces an expan¬ 
sion of the projected free energy that includes either a 
quadratic term or a term after the linear term that has 
the form |p|^, with C > 1- 


9.2. Attractions at Periodic Ranges of Separation 
via Oscillatory Priedel-Type Interactions 


Oscillatory (in sign) interactions between steps, medi¬ 
ated by substrate conduction electrons, ipso facto lead 
to attractive interactions between steps. As reviewed 
by Einstein [1996], such interactions have been know for 
many decades to account for the ordered patterns of ad¬ 
sorbates on metal surfaces [Einstein and Schrieffer 1973]. 
While at short range all electrons contribute, asymptot¬ 
ically the interaction is dominated by the electron [s] at 
the Fermi surface or, from another perspective, the non- 
analyticity in the response function at the nesting vector. 
The interaction energy has the form 

^par'" cx cos{2kF£ + «>) (44) 


FIG. 25: Perspective views of essentially the ECS (actually 
the Andreev surface free energy divided by fcsT) around the 
(001) facet calculated by the transfer matrix method with 
the PWFRG algorithm at ksT/ei = 0.3. (a) p-RSOS model 
{eintjei — —0.5).. (b) For comparison, the original unsticky 
RSOS model {tint = 0) From Akutsu [2011] . 


This, or its analogue for interacting local magnetic 
species, is called the RKKY [Ruderman and Kittel 1954, 
Yosida 1957] interaction. (The community studying mag¬ 
netism now labels as RKKY any interaction mediated by 
substrate electrons, not just the asymptotic limit writ¬ 
ten down in the RKKY papers.) The phase factor $ is 
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the nonperturbative result is the scattering phase shifts 
at the two atoms that are interacting; it is absent in 
the perturbational approach to this problem used in the 
RKKY papers. The exponent n indicates the decay en¬ 
velope. For interacting bulk entities, n = 3, the standard 
RKKY results. On metal surfaces, the leading term in 
the propagator is canceled by the image charge, leading 
to n = 5, with very rapid decay [Einstein and Schrief- 
fer 1973, Einstein 1996]. Such effects are insignificant 
for adatom interactions but can be more potent when 
a whole step participates. Redfield and Zangwill [1992] 
show that a line of localized perturbations will generate 
an interaction with n reduced by subtracting 1/2 and 
$ augmented by 7r/4. They used this result, with n = 
9/2, to account for Frohn et al. [1991] remarkable exper¬ 
imental results on vicinal Cu(OOl): from their observed 
bimodal TWD, Frohn et al. deduced that the step-step 
interaction is attractive for intermediate distances 3-5 
atoms. Indeed, it was their striking observation that led 
to several of the previously-discussed theory papers that 
claimed to find long-range step attractions. 

When there are metallic surface states (i.e. surface 
states for which their 2D band dispersion relation crosses 
the Fermi energy Ep) of Shockley nature (lying in a 
2D band gap containing Ep) , the indirect interaction 
has a much slower decay, with n = 2 [Lau and Kohn 
1978, Einstein 1996, Repp et al. 2000, Hyldgaard and 
Persson 2000, Knorr et al. 2002, Hyldgaard and Ein¬ 
stein 2005] Eurthermore, the Eermi wavevector typically 
is much smaller than that of bulk states, so the period of 
the oscillation in real space is much larger. Perhaps the 
most familiar metallic surface on metals is that at the 
center of the surface Brillouin zone (f) of the (111) sur¬ 
faces of noble metals, which exist inside the necks of the 
Eermi surface, discussed in textbooks like Ashcroft and 
Mermin [1976]. This is the state produces the famous 
wave structure in Eigler’s group’s dramatic STM images 
[Crommie et al. 1993] of atoms on metal surfaces. How¬ 
ever, there is a less well known metallic surface state on 
Cu (001), discovered relatively late (compared to other 
surface states) by Kevan [1983], centered at the zone-edge 
center X rather than F, that may provide a better expla¬ 
nation of Frohn et al.’s results in the Redfield-Zangwill 
framework. For surface-state mediated interactions be¬ 
tween steps, their formula indicates n = 3/2, comparable 
to the entropic and elastic repulsions. 

The effect of surface-state mediated interactions on 
TWDs was elucidated by Pai et al. [1994] in com¬ 
bined experimental and theoretical examination of vic¬ 
inal Ag(llO), which has a metallic surface state centered 
at Y, the middle of the shorter edge of the rectangular 
surface Brillouin zone [Liu et al. 1984]. In essence, the 
surface state introduces a second length scale, the Fermi 
wavelength Xp, in addition to (£), with the major con¬ 
sequence that the TWD is no longer a function of the 
single scaled dimensionless variable s but depends also 
on {£). With a suitable model potential Pai et al. [1994] 
could then account for the different TWDs at a few differ¬ 


ent misorientations (i.e., mean step spacings). Indeed, to 
establish convincingly that this Friedel-like effect is sig¬ 
nificant, one must measure several different values of {£). 
While this paper has been cited with regards to other 
modifications of TWDs (cf. e.g., Mugarza et al. [2006] 
and Li et al. [2010]), I have found no other investigations 
of Friedel-like effects on TWDs for several misorienta¬ 
tions of the same substance. 

Patrone and Einstein [2012] discuss other issues related 
to possible anisotropic surface state dispersion as well 
as showing the insensitivity to the point in the surface 
Brillouin zone about which the state is centered. 


10. CONCLUSIONS 

An aspect of ECS studies on which there has been 
substantial progress since the 1980’s, but which has re¬ 
ceived little attention in this chapter, is comparing and 
reconciling the values of the characteristic energies (sur¬ 
face free energy per area, step free energy per length, 
and step-step repulsion strength) that are extracted from 
experimental measurements with ever-improving calcu¬ 
lations (using density functional theory) of these ener¬ 
gies. Bonzel’s review [2003],as well as Nowicki and Bonzel 
[2004], Yu et al. [2006], Bonzel et al. [2007],Barreteau et 
al. [2003] contain extensive coverage of this issue for the 
soft metals to which his group has devoted exhaustive at¬ 
tention. Jeong and Williams [1999] review most results 
for silicon. Such efforts to find absolute energies has also 
taken place in studies of island shapes, e.g. of TiN (001) 
[Kodambaka et al. 2002] and (111) [Kodambaka et al. 
2003]. 

There are several significant advances in generic un¬ 
derstanding of ECS since the 1980’s. The Pokrovsky- 
Talapov {£) = 3/2) critical phenomena near the edge of 
the smoothly-curved region near a facet has proved to 
be far more robust and general than originally realized, 
while novel behavior is predicted in a very special direc¬ 
tion. Even though invoked in many accounts of sharp 
edges, long-range attractive £~'^ do not have an appar¬ 
ent physical basis, except perhaps in idiosyncratic cases. 
The likely cause is a reconstruction or adsorption that 
changes the surface free energy of the facet orientation. 
On the other hand, hill and valley structures are widely 
seen, and the possibility of azimuthal in addition to polar 
misorientation can lead to astonishingly rich phase dia¬ 
grams. Of course, non equilibrium considerations open 
up a whole new universe of behavior. Furthermore, at the 
nanoscale cluster shape is very sensitive to the particulars 
of a system, with the addition or removal of a single atom 
leading to a substantial change in shape, rather like bio¬ 
logical systems, in contrast to the macroscale phenomena 
that have been treated in this chapter. 



24 


Acknowledgements 

My research related to this subject was long supported 
by the UMD-NSF MRSEC; it is now supported par¬ 
tially by NSF-CHE 13-05892. Much of this paper is 
based on extensive collaboration with the experimental 
surface physics group at U. of Maryland, led by Ellen 
D. Williams until 2010, with ongoing guidance by Jan¬ 
ice Reutt-Robey and William G. Cullen, as well as with 
Margret Giesen and Harald Ibach, and Hans P. Bonzel, 
at FZ-Jtilich, partially sponsored by a Humboldt Senior 
U.S. Scientist Award, and, over the last decade, with 
Alberto Pimpinelli. 1 also benefited from theory inter¬ 
actions at Maryland with John D. Weeks and Dioni- 
sios Margetis, postdocs Olivier Pierre-Louis, Howard L. 
Richards, Ferenc Szalma, and Kwangmoo Kim, and grad¬ 
uate students Norman C. Bartelt, Raymond C. Nelson, 
Sanjay V. Khare, Hailu Gebremarian, Timothy J. Stase- 
vich, Rajesh Sathiyanarayan, Paul N. Patrone, and Josue 
R. Morales-Cifuentes. 


Bibliography 

D. B. Abraham and P. Reed, Phase Separation in the 
Two-Dimensional Ising Ferromagnet, Phys. Rev. Lett. 
33 (1974) 377. 

D. B. Abraham and P. Reed, Interface profile of the Ising 
ferromagnet in two dimensions, Commun. Math. Phys. 
49 (1976) 35. 

Y. Akutsu, N. Akutsu, and T. Yamamoto, Universal 
Jump of Gaussian Curvature at the Facet Edge of a 
Crystal, Phys. Rev. Lett. 61 (1988) 424. 

Y. Akutsu, N. Akutsu, and T. Yamamoto, Universality 
of the Tangential Shape Exponent at the Facet Edge 
of a Crystal, arXiv:cond-mat/9803189 (1998). 

Y. Akutsu and N. Akutsu, The Equilibrium Facet Shape 
of the Staggered Body-Centered-Cubic Solid-on-Solid 
Model, Prog. Theor. Phys. 116 (2006) 983. 

N. Akutsu, Non-universal equilibrium crystal shape re¬ 
sults from sticky steps, J. Phys.: Condens. Matter 23 
(2011) 485004. 

O. L. Alerhand, D. Vanderbilt, R. D. Meade, and J. D. 
Joannopoulos, Spontaneous Formation of Stress Do¬ 
mains on Crystal Surfaces, Phys. Rev. Lett. 61 (1988) 
1973. 

F. Almgren and J. E. Taylor, Optimal Geometry in Equi¬ 
librium and Growth, Fractals 3 (1996) 713. 

A. F. Andreev, Faceting phase transitions of crystals, 
Sov. Phys.-JETP 53 (1982) 1063 [Zh. Eksp. Teor. Phys 
80 (1981) 2042. 

N.W. Ashcroft and N. D. Mermin, Solid State Physics 
(Cengage Learning, 1976). 

J. E. Avron, H. van Beijeren, L. S. Schulman, and R. K. P. 
Zia, Roughening transition, surface tension and equi¬ 
librium droplet shapes in a two-dimensional Ising sys¬ 
tem, J. Phys. A: Math. Gen. 15 (1982) L81. 

S. Balibar and B. Castaing, J. Phys. Lett. 41 (1980) L329. 


S. Balibar, H. Alles, and A. Ya. Parshin, The surface of 
helium crystals. Rev. Mod. Phys. 77 (2005) 317. 

C. Barreteau, F. Raouafi, M. C. Desjonqueres, and D. 
Spanjaard, J. Phys.: Condens. Matter 15 (2003) 3171. 

N. C. Bartelt, E. D. Williams, R. J. Phaneuf, Y. Yang, 
and S. Das Sarma, Orientational stability of silicon 
surfaces, J. Vac. Sci. Technol. A 7 (1989) 1898. 

N. C. Bartelt, T. L. Einstein, and E. D. Williams, The 
Influence of Step-Step Interactions on Step Wandering, 
Surface Sci. 240 (1990) L591. 

A. A. Baski, S. C. Erwin, and L. J. Whitman, The struc¬ 
ture of silicon surfaces from (001) to (111), Surface Sci. 
392 (1997) 69. 

J. M. Bermond, J. J. Metois, J. C. Heyraud, and F. Flo¬ 
ret, Shape universality of equilibrated silicon crystals. 
Surface Science 416 (1998) 430. 

F. Bonczek, T. Engel, and E. Bauer, Surface Sci. 97 
(1980) 595. 

H. P. Bonzel and E. Preuss, Morphology of periodic sur¬ 
face profiles below the roughening temperature: as¬ 
pects of continuum theory, Surface Sci. 336 (1995) 209. 

H. P. Bonzel, 3D equilibrium crystal shapes in the new 
light of STM and AFM, Physics Reports 385 (2003) 1. 

H. P. Bonzel, D.K. Yu, and M. Schefher, The three- 
dimensional equilibrium crystal shape of Pb: Recent 
results of theory and experiment, Appl. Phys. A 87 
(2007) 391. 

H. P. Bonzel and M. Nowicki, Absolute surface free en¬ 
ergies of perfect low-index orientations of metals and 
semiconductors, Phys. Rev. B 70 (2004) 245430. 

W. K. Burton, N. Gabrera, and F. C. Frank, The growth 
of crystals and the equilibrium structure of their sur¬ 
faces, Phil. Trans. Roy. Soc. A 243 (1951) 299. 

N. Gabrera and N. Garci'a, Roughening transition in the 
interface between superfluid and solid "^He, Phys. Rev. 
B 25 (1982) 6057. 

N. Gabrera, The equilibrium of crystal surfaces. Surface 
Sci. 2 (1964) 320. 

J.W. Cahn, and W. C. Carter, Crystal Shapes and 
Phase Equilibria: A Common Mathematical Basis, 
Metallurgical and Materials Transactions A-Physical 
Metallurgy and Materials Science 27 (1996) 1431. 
arXiv:cond-mat/0703564vl] . 

J.W. Cahn, Transitions and Phase-Equilibria Among 
Grain-Boundary Structures, J. de Physique, 43(C6) 
(1982) 199, Proceedings of Conference on the Struc¬ 
ture of Grain Boundaries, Caen, France. 

H.B. Callen, Thermodynamics and an Introduction to 
Themostatistics, 2nd ed. (Wiley, New York, 1985). 

F. Calogero, Solution of a Three-Body Problem in One 
Dimension, J. Math. Phys. 10 (1969) 2191, Ground 
State of a One-Dimensional A-Body System, J. Math. 
Phys. 10 (1969) 2197. 

E. Carlon and H. van Beijeren, Equilibrium shapes and 
faceting for ionic crystals of body-centered-cubic type, 
Phys. Rev. E 62 (2000) 7646. 

Y. Carmi, S. G. Lipson, and E. Polturak, Critical behav¬ 
ior of vicinal surfaces of "^He, Phys. Rev. B 36 (1987) 



25 


1894. 

R. Cerf and J. Picard, The Wulff Crystal in Ising and 
Percolation Models (Springer, Berlin, 2006). 

W.-C. Cheng and P. Wynblatt, Coupled compositional 
and roughening phase transitions at the surface of a 
Pb-Bi-Ni alloy. Surface Sci. 302 (1994) 185. 

A. A. Chernov, Layer-Spiral Growth of Crystals, Sov. 
Phys.-Uspekhi 4 (1961) 116. 

C.V. Ciobanu, D.T. Tambe, V. B. Shenoy, C. Z. Wang, 
and K. M. Ho, Atomic-scale perspective on the origin 
of attractive step interactions on Si(113), Phys. Rev. 
B 68 (2003) 201302R. 

M. Copel and T. Gustafsson, Structure of Au(llO) De¬ 
termined with Medium-Energy-Ion scattering, Phys. 
Rev. Lett. 57 (1986) 723. 

M. F. Crommie, C. P. Lutz, and D. M. Eigler. Imaging 
standing waves in a two-dimensional electron gas, Na¬ 
ture 363 (1993) 524; Confinement of electrons to quan¬ 
tum corrals on a metal surface. Science 262 (1993) 218. 

P. Curie, Bull. Soc. Min. de France 8 (1885) 145. 

B. Dacorogna and C.E. Pfister, Wulff theorem and best 
constant in Sobolev inequality, J. Math. Pures. Appl. 
71 (1992) 97. 

S. R. Dahmen, B. Wehefritz, and G. Albertini, A 
novel exponent in the Equilibrium Shape of Grystals, 
arXiv:cond-mat/9802152 (1998). 

M. S.Daw, S. M.Foiles, and M. I.Baskes, The Embedded- 
Atom Method - a Review of Theory and Applications, 
Mater. Sci. Rep. 9 (1993) 251. 

J. De Coninck, F. Dunlop, and V. Rivasseau, On the mi¬ 
croscopic validity of the Wulff construction and of the 
generalized Young equation, Commun. Math. Phys. 
121 (1989) 401. 

P. G. de Gennes, Soluble Model for Fibrous Structures 
with Steric Constraints, J. Chem. Phys. 48 (1968) 
2257. 

M. Degawa, T. J. Stasevich, W. G. Cullen, A. Pimpinelli, 

T.L. Einstein, and E. D. Williams, Distinctive Fluc¬ 
tuations in a Confined Geometry, Phys. Rev. Lett. 97 
(2006) 080601. 

M. Degawa, T. J. Stasevich, A. Pimpinelli, T. L. Einstein, 
and E. D.Williams, Facet-edge Fluctuations with Pe¬ 
riphery Diffusion Kinetics, Surface Sci. 601 (2007) 3979 
[Proc. EGOSS 2006]. 

For a review of fermionic methods, see M. den 
Nijs, “The Domain Wall Theory of Two-dimensional 
Commensurate-Incommensurate Phase Transitions,” 
in: Phase Transitions and Critical Phenomena, Vol. 
12, edited by C. Domb and J. L. Lebowitz (Academic, 
London, 1989) pp. 219-333. 

S. Dieluweit, H. Ibach, M. Giesen, and T.L. Einstein, 
Orientation Dependence of Step Stiffness: Failure of 
Solid-on-Solid and Ising Models to Describe Experi¬ 
mental Data, Phys. Rev. B 67 (2003) 121410(R). 

S. Dietrich and A. Haase, Scattering of X-rays and neu¬ 
trons at interfaces, Physics Reports 260 (1995) 1. 

A. Dinghas, Uber einen Gcometrischen Satz von Wulff fiir 
die Gleichgewichtsform von Kristallen, Z. Kristallog. 


105 (1944) 304. 

R. L. Dobrushin, R. Kotecky, and S. B. Shlosman, A Mi¬ 
croscopic Justification of the Wulff Construction, J. 
Stat. Phys. 12 (1993) 1. 

R. L. Dobrushin, R. Kotecky, and S. Shlosman, 
The Wnlff Construction: A Global Shape from 
Local Interactions (AMS, Providence, Rhode 
Island, 1992); preprint: http://www.cpt.univ- 

mrs.fr/dobrushin/DKS-book.pdf, 

P. M. Duxbury and T. J. Pence (eds.). Dynamics of Crys¬ 
tal Surfaces and Interfaces (Plenum, New York, 1997). 

F. J. Dyson, Statistical Theory of the Energy Levels of 
Complex Systems. HI, J. Math. Phys. 3 (1962) 166. 

F. J. Dyson, Correlations between eigenvalues of a ran¬ 
dom matrix, Commun. Math. Phys. 19 (1970) 235. 

T. L. Einstein and J. R. Schrieffer, Indirect Interaction 
Between Adatom Pairs on a Tight-Binding Solid, 
Phys. Rev. B7 (1973) 3629. 

T. L. Einstein, “Interactions Between Adsorbate Parti¬ 
cles,” in: Physical Structure of Solid Surfaces, edited 
by W. N. Unertl (Elsevier, Amsterdam, 1996), Hand¬ 
book of Surface Science, vol. 1, S. Holloway and N. V. 
Richardson, series eds., pp. 577-650. 

T. L. Einstein, H. L.Richards, S. D.Cohen, and O. Pierre- 
Louis, Terrace-Width Distributions and Step-Step Re¬ 
pulsions on Vicinal Surfaces: Symmetries, Scaling, 
Simplifications, Subtleties, and Schrodinger, Surface 
Sci. 493 (2001) 460. 

T. L. Einstein, Using the Wigner-Ibach Surmise to An¬ 
alyze Terrace-Width Distributions: History, User’s 
Guide, and Advances, Appl. Phys. A 87 (2007) 375 
cond-mat/0612311|. 

T. L. Einstein and Alberto Pimpinelli, Dynamical Scal¬ 
ing Implications of Ferrari, Prhofer, and Spohn’s Re¬ 
markable Spatial Scaling Results for Facet-Edge Fluc¬ 
tuations, J. Statistical Phys. 155 (2014) 1178; longer 
version posted at arXiv 1312.4910vl. 

A. Emundts, H.P. Bonzel, P. Wynblatt, K. Thrmer, 
J. Reutt-Robey, and E. D. Williams, Gontinuous and 
discontinuous transitions on 3D equilibrium crystal 
shapes: a new look at Pb and Au, Surface Sci. 481 
(2001) 13. 

F. Ercolessi, A. Bartolini, M. Garofalo, M. Parrinello, 
and E. Tosatti, Au Surface Reconstructions in the Glue 
Model, Surface Sci 189/190 (1987) 636. 

P. L. Ferrari and H. Spohn, Step Fluctuations for a 
Faceted Crystal, J. Stat. Phys. 113 (2003) 1. 

P. L. Ferrari, M. Prahofer, and H. Spohn, Fluctuations 
of an atomic ledge bordering a crystalline facet, Phys. 
Rev. E 69 (2004) 035102. 

D. S. Fisher and J. D. Weeks, Shape of Crystals at 
Low Temperatures: Absence of Quantum Roughening, 
Phys. Rev. Lett. 50 (1983) 1077. 

M. P. A. Fisher, D. S. Fisher, and J. D. Weeks, Agreement 
of Capillary-Wave Theory with Exact Results for the 
Interface Profile of the Two-dimensional Ising Model, 
Phys. Rev. Lett. 48 (1982) 368. 


26 


M. Flytzani-Stephanopoulos and L. D. Schmidt, Mor¬ 
phology and Etching Processes on Macroscopic Metal 
Catalysts, Prog. Surf. Sci. 9 (1979) 83. 

I. Fonseca and S. Muller, An uniqueness proof for the 
Wulff problem, Proc. Royal Soc. Edinburgh 119 (1991) 
125. 

I. Fonseca, The Wulff theorem revisited, Proc. R. Soc. 
London A 432 (1991) 125. 

F. C. Frank, Metal Surfaces, (American Society for Met¬ 
als, Metals Park, Ohio 1962/3) chap. 1, pp. 1-15. 

J. Frohn, M. Giesen, M. Poensgen, J. F. Wolf, and 

H. Ibach, Attractive interaction between steps, Phys. 
Rev. Lett. 67 (1991) 3543. 

G. Gallavotti, The phase separation line in the two- 
dimensional Ising model, Commun. Math. Phys. 27 
(1972) 103. 

N. Garcia, J. J. Saenz, and N. Cabrera, Cusp Points in 
Surface Free Energy: Faceting and First Order Phase 
Transitions, Physica 124B (1984) 251. 

Hailu Gebremariam, S. D. Cohen, H. L. Richards, and 
T.L. Einstein, Analysis of Terrace Width Distribu¬ 
tions Using the Generalized Wigner Surmise: Galibra- 
tion Using Monte Garlo and Transfer-Matrix Galcula- 
tions, Phys. Rev. B 69 (2004) 125404. 

J. W. Gibbs, On the Equilibrium of Heterogeneous Sub¬ 
stances, Transactions of the Gonnecticut Academy of 
Arts and Sciences, 3, 108248, 343524, (18741878). Re¬ 
produced in both The Scientific Papers (1906), pp. 
55-353 and The Gollected Works of J. Willard Gibbs, 
edited by W. R. Longley and R. G. Van Name (Yale 
University Press, New Haven, 1957 [1928]) pp. 55-353. 

J.W. Gibbs, Trans Gonnecticut Acad. 3, 108-248 (1877) 
343-524; see Collected Works 1928 (previous entry). 

M. Giesen, Step and island dynamics at solid/vacuum 
and solid/liquid interfaces. Prog. Surface Sci. 68 (2001) 

I . 

J. Gladic, Z. Vucic, and D. Lovric, Gritical behaviour of 
the curved region near 111-facet edge of equilibrium 
shape cuprous selenide large single crystals, J. Grystal 
Growth 242 (2002) 517532. 

E. E. Gruber and W.W. Mullins, J. Phys. Ghem. Solids 
28 (1967) 875. 

T. Guhr, A. Miiller-Groeling, and H. A. Weidenmiiller, 
Random-matrix theories in quantum physics: common 
concepts, Phys. Rept 299 (1998) 189. 

F. Haake, Quantum Signatures of Ghaos, 2nd ed. 
(Springer, Berlin, 1991). 

F. D. M. Haldane and J. Villain, J. Phys. (Paris) 42 
(1981) 1673. 

C.A. Handwerker, M. D. Vaudin, and J.E. Blendell, 
Equilibrium Crystal Shapes and Surface Phase Dia¬ 
grams at Surfaces in Ceramics, J. Phys. Colloques 49 
(1988) C5-367. 

P. Hecquet, Surface stress modifies stepstep interaction 
energy with respect to the MarchenkoParshin model, 
Surface Sci. 602 (2008) 819. 

G. A.Held, D. M.Goodstein, and J. D.Brock, Phase sepa¬ 
ration and step roughening of vicinal Si(ill): An x-ray¬ 


scattering study, Phys. Rev. B 51 (1995) 7269. 

C. Herring, Some Theorems on the Free Energies of Crys¬ 
tal Surfaces, Phys. Rev. 82 (1951) 87. 

C. Herring, The use of classical macroscopic concepts in 
surface energy problems, in: Structure and Properties 
of Solid Surfaces, edited by R. Corner and C. S.Smith 
(University of Chicago Press, Chicago, 1953), Chapter 

I, pp. 5-81. 

H. Hilton, Mathematical Crystallography (Oxford, 1903). 

D. W. Hoffman and J.W. Cahn, A Vector Thermody¬ 
namics for Anisotropic Surfaces: 1. Fundamentals and 
Application to Plane Surface Junctions, Surface Sci. 
31 (1972) 368; J.W. Cahn and D.W. Hoffman, Vector 
Thermodynamics for Anisotropic Surfaces: 11. Curved 
and Faceted Surfaces, Acta Met. 22 (1974) 1205. 

R. Q. Hwang, E. D. Williams, and R. L. Park, High- 
resolution low-energy electron-diffraction study of the 
phase diagram of vicinal Si(lll) surfaces, Phys. Rev. 
B 40 (1989) 11716. 

P. Hyldgaard and T. L. Einstein, Interactions Mediated 
by Surface States: From Pairs and Trios to Adchains 
and Ordered Overlayers, J. Crystal Growth 275 (2005) 
el637 cond-mat/0408645 . 

P. Hyldgaard and M. Persson, Long-ranged adsorbatead- 
sorbate interactions mediated by a surface-state band, 

J. Phys.: Condes. Matter 12 (2000) L13. 

H. Ibach and W. Schmickler, Step Line Tension on a 
Metal Electrode, Phys. Rev. Lett. 91 (2003) 016106. 

T. Ihle, G. Misbah, and O. Pierre-Louis, Equilibrium step 
dynamics on vicinal surfaces revisited, Phys. Rev. B 58 
(1998) 2289. 

K. A. Jackson, “Theory of Melt Growth,” in: Crys¬ 
tal Growth and Characterization, edited by R. Ueda 
and J. B. Mullin (North-Holland, Amsterdam, 1975) 
[Proc. lSSCG-2 Spring School, Lake Kawaguchi, Japan 
(1974)], pp. 21-32. 

C. Jayaprakash and W. F. Saam, Thermal evolution of 
crystal shapes: The fee crystal, Phys. Rev. B 30 (1984) 
3916. 

C. Jayaprakash, C. Rottman, and W. F. Saam, Simple 
model for crystal shapes: Step-step interactions and 
facet edges, Phys. Rev. B 30 (1984) 6549; in the Hamil¬ 
tonian in their Eq. 3, the factor t/2 should have been 
t. See Williams et al. [1994]. 

C. Jayaprakash, W. F. Saam, and S. Teitel, Roughening 
and Facet Formation in Crystals, Phys. Rev. Lett. 50 
(1983) 2017. 

H.-C. Jeong and E. D. Williams, Steps on Surfaces: Ex¬ 
periment and Theory, Surface Sci. Rept. 34 (1999) 171. 

H.-C. Jeong and J. D. Weeks, Effects of stepstep inter¬ 
actions on the fluctuations of an individual step on a 
vicinal surface and its wavelength dependence, Surface 
Sci. 432 (1999) 101, and references therein. 

B. Jobs, T.L. Einstein, and N. C.Bartelt, Distribution 
of Terrace Widths on a Vicinal Surface in the One- 
Dimensional Free-Fermion Model, Phys. Rev. B 43 
(1991) 8153. 



27 


M. Kalff, G. Comsa, and T. Michely, How Sensitive is 
Epitaxial Growth to Adsorbates?, Phys. Rev. Lett. 81 
(1998) 1255. 

M. Kardar and D. R. Nelson, Gommensurate- 
Incommensurate Transitions with Quenched Random 
Impurities, Phys. Rev. Lett. 55 (1985) 1157. 

K. O. Keshishev, A. Ya. Parshin, and A. V. Babkin, Grys- 
tallization Waves in He-4, Sov. Phys. JETP 53 (1981) 
362. 

S. D. Kevan, Observation of a new surface state on 
Gu(OOl), Phys. Rev. B 28 (1983) 2268(R). 

N. Knorr, H. Brune, M. Epple, A. Hirstein, M. A. Schnei¬ 
der, and K. Kern, Long-range adsorbate interactions 
mediated by a two-dimensional electron gas, Phys. 
Rev. B 65 (2002) 115420. 

S. Kodambaka, S.V. Khare, V. Petrova, A. Vailionis, 

I. Petrov, and J. E. Greene, Absolute orientation- 
dependent TiN(OOl) step energies from two- 
dimensional equilibrium island shape and coarsening 
measurements on epitaxial TiN(OOl) layers. Surface 
Sci. 513 (2002) 468. 

S. Kodambaka, S.V. Khare, V. Petrova, D.D. John¬ 
son, 1. Petrov, and J. E. Greene,. Absolute orientation- 
dependent anisotropic TiN(lll) island step energies 
and stiffnesses from shape fluctuation analyses, Phys. 
Rev. B 67 (2003) 035409. 

S. Kodambaka, S.V. Khare, 1. Petrov, J.E. Greene, 
Two-dimensional island dynamics: Role of step energy 
anisotropy, Surface Sci. Reports 60 (2006) 55. 

W. Kossel, Extending the Law of Bravais, Nachr. Ges. 
Wiss. Gottingen (1927) 143. 

W. Kossel, Zur Energetik von Oberflachenvorgangen, An¬ 
nul. Physik 21 (1934) 457. 

J. M. Kosterlitz and D. J. Thouless, Ordering, metastabil¬ 
ity and phase transitions in two-dimensional systems, 

J. Phys. C6 (1973) 1181. 

J. M. Kosterlitz, The critical properties of the two- 
dimensional xy model, J. Phys. C 7 (1974) 1046. 

R. V. Kukta, A. Peralta, and D. Kouris, Elastic Interac¬ 
tion of Surface Steps: Effect of Atomic-Scale Rough¬ 
ness, Phys. Rev. Lett. 88 (2002) 186102. 

M. Lassig, Vicinal Surfaces and the Calogero-Sutherland 
Model, Phys. Rev. Lett 77 (1996) 526. 

L. D. Landau and E. M. Lifshitz, Statistical Physics, Part 
1, 3rd edition revised and enlarged by E. M. Lifshitz 
and L. P. Pitaevskii (Pergamon Press, Oxford, 1980), 
§155. 

K. H. Lau and W. Kohn, Indirect Long-Range Oscillatory 
Interaction Between Adsorbed Atoms, Surface Sci. 75 
(1978) 69; T. L. Einstein, Gomment on K. H.Lau and 
W. Kohn: “Oscillatory Indirect Interaction between 
Adsorbed Atoms”-Non-Asymptotic Behavior in Tight- 
Binding Models at Realistic Parameters, Surface Sci. 
75 (1978) L161. 

H. J. Leamy, G.H. Gilmer, and K.A. Jackson, “Statis¬ 
tical Thermodynamics of Glean Surfaces,” in: Surface 
Physics of Materials, edited by J. M. Blakely, vol. 1 


(Academic Press New York, 1975), chap. 3, pp. 121- 
188. 

F. Li, F. Allegretti, S. Surnev, and F.P. Netzer, Atomic 
engineering of oxide nanostructure superlattices. Sur¬ 
face Sci. 604 (2010) L43. 

E. H. Lieb, Exact Solution of the Two-Dimensional Slater 
KDP Model of a Ferroelectric, Phys. Rev. Lett. 19 
(1967) 108. 

E. H. Lieb and F. Y. Wu, in: Phase Transitions and Crit¬ 
ical Phenomena, edited by C. Domb and M. S. Green 
(Academic Press, London, 1972) vol. 1, p. 331. 

H. Liebmann, Der Curie-Wulff’sche Satz iiber Combina- 
tionsformen von Krystallen, Z. Kristallog. 53 (1914) 
171. 

S. H. Liu, C. Hinnen, C. N. van Huong, N. R. de Tacconi, 
and K.-M. Ho, Surface State Effects on the Electrore¬ 
flectance Spectroscopy of Au Single Crystal Surfaces, 
J. Electroanal. Chem. 176 (1984) 325. 

F. Liu, “Modeling and Simulation of Strain-mediated 
Nanostructure Formation on Surface,” in: Handbook 
of Theoretical and Computational Nanotechnology, 
vol. 4, edited by Michael Rieth and Wolfram Schom- 
mers, American Scientific Publishers (2006), pp. 577- 
625. 

I. Lyuksyutov, A. G. Naumovets, and V. Pokrovsky, Two- 
Dimensional Crystals (Academic Press, San Diego, 
1992). 

T. E. Madey, C.-H. Nien, K. Pelhos, J. J. Kolodziej, 1. M. 
Abdelrehim, and H.-S. Tao, Faceting induced by ultra- 
thin metal films: structure, electronic. Surface Sci. 438 
(1999) 191. 

T.E. Madey, J. Guan, C.-H. Nien, H.-S. Tao, C.-Z. Dong, 
and R. A. Campbell, Surf. Rev. Lett. 3 (1996) 1315. 

V. 1. Marchenko and A. Ya. Parshin, Elastic properties of 
crystal surfaces, Sov. Phys. JETP 52 (1980) 129 [Zh. 
Eksp. Teor. Fiz. 79 (1980) 257.] 

J. J. Metois and J. C. Heyraud, Analysis of the critical 
behaviour of curved regions in equilibrium shapes of 
In crystals. Surface Sci. 180 (1987) 647. 

M. L. Mehta, Random Matrices, 3rd ed. (Academic, New 
York, 2004). 

H. Meltzman, D. Chatain, D. Avizemer, T. M. Besmann, 
and W. D. Kaplan, The equilibrium crystal shape of 
nickel, Acta Materialia 59 (2011) 3473. 

T. Michely and J. Krug, Islands, Mounds, and Atoms: 
Patterns and Process in Crystal Growth Far from 
Equilibrium, Springer, Berlin, 2004, chap. 3. 

S. Miracle-Sole, Facet Shapes in a Wulff Crystal, 
in: Mathematical Results in Statistical Mechanics, 
edited by S. Miracle-Sole, J. Ruiz, and V. Zagrebnov 
(World Scientific, Singapore, 1999) pp. 83-101 [arXiv: 
1206.3736vl]. 

S. Miracle-Sole and J. Ruiz, On the Wulff construction 
as a problem of equivalence of statistical ensembles, 
in: On Three Levels: Micro, Meso and Macroscopic 
Approaches in Physics, edited by M. Fannes and A. 
Verbeure (Plenum Press, New York, 1994), pp. 295- 
302 arXiv:1206.3739vl]. 


28 


S. Miracle-Sole, WulfF shape of equilibrium crystals 
(2013) [arXiv: 1307.5180vl]. 

A. Mugarza, F. Schiller, J. Kuntze, J. Cordon, M. Ruiz- 
Oses, and J. E. Ortega, Modelling nanostructures with 
vicinal surfaces, J. Phys.: Condens. Matter 18 (2006) 
S27. 

B. Muller, L. Nedelmann, B. Fischer, H. Brune, J.V. 
Barth, and K. Kern, Island Shape Transition in Het- 
eroepitaxial Metal Growth on Square Lattices, Phys. 
Rev. Lett. 80 (1998) 2642. 

W. W. Mullins, “Solid Surface Morphologies Governed by 
Capillarity,” in: Metal Surfaces: Structure, Energet¬ 
ics, and Kinetics, ed. by W. D. Robertson and N. A. 
Gjostein (American Society for Metals, Metals Park 
(OH), 1962/3), chap. 2, pp. 17-62. 

R. Najafabadi and D. J. Srolovitz, Elastic step interac¬ 
tions on vicinal surfaces of fee metals. Surface Sci. 317 
(1994) 221. 

R. C. Nelson, T. L. Einstein, S. V. Khare, and P. J. Rous, 
Energies of Steps, Kinks, and Defects on AglOO and 
111 Using Embedded Atom Method, and Some Con¬ 
sequences, Surface Sci. 295 (1993) 462. 

T. Nishino and K. Okunishi, J. Phys. Soc. Jpn. 64 (1995) 
4084. 

D.Y. Noh, K. 1. Blum, M. J. Ramstad, and R. J. Bir- 
geneau, Long-range coherence and macroscopic phase 
separation of steps on vicinal Si(lll), Phys. Rev. B 44 
(1991) 10969. 

D.Y. Noh, K. 1. Blum, M. J. Ramstad, and R. J. Birge- 
neau. Faceting, roughness, and step disordering of vic¬ 
inal Si(lll) surfaces: An x-ray-scattering study, Phys. 
Rev. B 48 (1993) 1612. 

D.Y. Noh, K. S. Liang, Y. Hwu, and S. Chandavarkar, 
Surface Sci. 326 (1995) L455. 

1. M. Nolden and H. van Beijeren, Equilibrium shape of 
bcc crystals: Thermal evolution of the facets, Phys. 
Rev. B 49 (1994) 17224. 

M. Nowicki, C. Bombis, A. Emundts, H. P. Bonzel, and 
P. Wynblatt, Universal exponents and step-step inter¬ 
actions on vicinal Pb(lll) surfaces, Eur. Phys. Lett. 
59 (2002a) 239. 

M. Nowicki, C. Bombis, A. Emundts, H. P. Bonzel, and 
P. Wynblatt, Stepstep interactions and universal expo¬ 
nents studied via three-dimensional equilibrium crys¬ 
tal shapes. New J. Phys. 4 (2002b) 60. 

M. Nowicki, C. Bombis, A. Emundts, and H. P. Bonzel, 
Absolute step and kink formation energies of Pb de¬ 
rived from step roughening of two-dimensional islands 
and facets, Phys. Rev. B 67 (2003) 075405. 

P. Nozieres, Shape and Growth of Crystals, in: Solids Far 
From Equilibrium, edited by C. Godreche, Cambridge 

U. Press, Cambridge, 1992. 

B. Ocko, private discussion, 2014. 

K. Okunishi, Y. Hieida and Y. Akutsu, i5-function Bose- 
gas picture of S=1 antiferromagnetic quantum spin 
chains near critical fields, Phys. Rev. B 59 (1999) 6806; 
Middle-field cusp singularities in the magnetization 


process of one-dimensional quantum antiferromagnets, 
ibid. 60 (1999) R6953. 

B. Olshanetsky and V. Mashanov, LEED Studies of 
Clean High Miller Index Surfaces of Silicon, Surface 
Sci. Ill (1981) 414. 

W. W. Pai, J. S. Ozcomert, N. C. Bartelt, T. L. Einstein, 
and J. E. Reutt-Robey, Terrace-Width Distributions 
on Vicinal Ag(llO): Evidence of Oscillatory Interac¬ 
tions, Surface Sci. 307-309 (1994) 747. 

P. N. Patrone and T. L. Einstein, Anisotropic Surface 
State Mediated RKKY Interaction Between Adatoms, 
Phys. Rev. B 85 (2012) 045429. 

R. F. Peierls, Quantum Theory of Solids (Clarendon, Ox¬ 
ford, 1955), p. 108. (Also H. Frohlich, Proc. Roy. Soc. 
A 223 (1954) 296. 

D. Peng, S. Osher, B. Merriman, and H.-K. Zhao, “The 
geometry of WulfF crystal shapes and its relation with 
Riemann problems,” in: Nonlinear Partial Differential 
Equations, Contemp. Math. 238, AMS, Providence, 
RI, 1999, pp. 251-303. 

C. E. Pfister, Large deviations and phase separation in 
the two-dimensional Ising model, Helv. Phys. Acta 64 
(1991) 953. 

R. J. Phaneuf and E. D. Williams, Surface Phase Sepa¬ 
ration of Vicinal Si (111), Phys. Rev. Lett. 58 (1987) 
2563. 

R. J. Phaneuf, N. C. Bartelt, E. D. Williams, W. Swiech, 
and E. Bauer, Crossover from Metastable to Unstable 
Facet Growth on Si(lll), Phys. Rev. Lett. 71 (1993) 
2284. 

A. Pimpinelli and J. Villain, Physics of Grystal Growth 
(Gambridge University Press, Cambridge, England, 
1998). 

A. Pimpinelli, M. Degawa, T. L. Einstein, and E. D. 
Williams, A Facet Is Not an Island: Step-Step Interac¬ 
tions and the Fluctuations of the Boundary of a Crys¬ 
tal Facet, Surface Sci. 598 (2005) L355. 

V. L. Pokrovsky and A. L. Talapov, Ground State, Spec¬ 
trum, and Phase Diagram of Two-Dimensional Incom¬ 
mensurate Grystals, Phys. Rev. Lett. 42 (1979) 65; 
Sov. Phys.-JETP 51 (1980) 134. 

V. L. Pokrovsky and A. L. Talapov, Theory of Incommen¬ 
surate Grystals, Soviet Scientific Reviews Supplement 
Series Physics vol. I, Harwood Academic Publishers, 
Chur, 1984, and references therein. 

G. Prevot and B. Groset, Revisiting Elastic Interactions 
between Steps on Vicinal Surfaces: The Buried Dipole 
Model, Phys. Rev. Lett. 92 (2004) 256104. 

A. C. Redfield and A. Zangwill, Attractive interactions 
between steps, Phys. Rev. B 46 (1992) 4289. 

J. Repp, F. Moresco, G. Meyer, K.-H. Rieder, P. 
Hyldgaard, and M. Persson, Substrate Mediated Long- 
Range Oscillatory Interaction between Adatoms: Gu 
/Cu(lll), Phys. Rev. Lett. 85 (2000) 2981. 

C. Rottman and M. Wortis, Statistical mechanics of equi¬ 
librium crystal shapes: Interfacial phase diagrams and 
phase transitions, Phys. Rept. 103 (1984a) 59. 



29 


C. Rottman and M. Wortis, Equilibrium crystal shapes 
for lattice models with nearest- and next-nearest- 
neighbor interactions, Phys. Rev. B 29 (1984b) 328. 

C. Rottman, M. Wortis, J. C. Heyraud, and J. J. Metois, 
Equilibrium Shapes of Small Lead Crystals: Obser¬ 
vation of Pokrovsky-Talapov Critical Behavior, Phys. 
Rev. Lett. 52 (1984) 1009. 

S. Rousset, J. M. Berroir, V. Repain, Y. Garreau, V. H. 
Etgens, J. Lecoeur, and R. Pinchaux, Thermal faceting 
behavior of Au(4,5,5), Surface Sci. 443 (1999) 265. 

M. A. Ruderman and C. Kittel, Indirect Exchange Cou¬ 
pling of Nuclear Magnetic Moments by Conduction 
Electrons, Phys. Rev. 96 (1954) 99. 

J. J. Saenz and N. Garci'a, Classical Critical Behaviour 
in Crystal Surfaces Near Smooth and Sharp Edges, 
Surface Sci. 155 (1985) 24. 

R. Sathiyanarayanan, A. BH. Hamouda, and T. L. Ein¬ 
stein, Terrace-width Distributions of Touching Steps: 
Modihcation of the Fermion, with Implications for 
Measuring Step-step Interactions, Phys. Rev. B 80 
(2009) 153415. 

R. Sato and Y. Akutsu, Curvature Jump at the Facet 
Edge of a Crystal for Arbitrary Surface Orientation, 
J. Phys. Soc. Jpn. 64 (1995) 3593. 

R. E. Schlier and H. E. Farnsworth, Structure and Ad¬ 
sorption Characteristics of Clean Surfaces of Germa¬ 
nium and Silicon, J. Ghem. Phys. 30 (1959) 917. 

H. J. Schulz, B. 1. Halperin, and G. L. Henley, Dislo¬ 
cation interaction in an adsorbate solid near the 
commensurate-incommensurate transition, Phys. Rev. 
B 26 (1982) 3797. 

R. F. Sekerka, Theory of Crystal Growth Morphology, in: 
Grystal Growth - From Fundamentals to Technology, 
edited by G. Miiller, J.-J. Metois, and P. Rudolph (El¬ 
sevier, Amsterdam, 2004), pp. 55-93. 

V. B. Shenoy and C.V. Ciobanu, Ab initio density func¬ 
tional studies of stepped TaG surfaces, Phys. Rev. B 
67 (2003) 081402. 

N. Shimoni, S. Ayal, and O. Millo, Step dynamics 
and terrace-width distribution on flame-annealed gold 
films: The effect of step-step interaction, Phys. Rev. B 
62 (2000) 13147. 

G.A. Somorjai and M. A. Van Hove, Adsorbate-Induced 
Restructuring of Surfaces, Prog. Surf. Sci. 30 (1989) 
201 . 

S. Song and S. G. J. Mochrie, Tricriticality in the orienta¬ 
tional phase diagram of stepped Si(113) surfaces, Phys. 
Rev. Lett. 73 (1994) 995. 

S. Song and S. G. J. Mochrie, Attractive step-step inter¬ 
actions, tricriticality, and faceting in the orientational 
phase diagram of silicon surfaces between [113] and 
[114], Phys. Rev. B 51 (1995) 10068. 

S. Song, M. Yoon, and S. G. J. Mochrie, Faceting, tricriti¬ 
cality, and attractive interactions between steps in the 
orientational phase diagram of silicon surfaces between 
[113] and [55 12], Surface Sci. 334 (1995) 153. 

T. J. Stasevich and T. L. Einstein, Analytic Formulas 
for the Orientation Dependence of Step Stiffness and 


Line Tension: Key Ingredients for Numerical Modeling, 
[SIAM-]Multiscale Model. Simul. 6 (2007) 90. 

T. J. Stasevich, C. Tao, W. G. Cullen, E. D. Williams, and 

T. L. Einstein, Impurity Decoration for Crystal Shape 
Control: Ceo on Ag(lll), Phys. Rev. Lett. 102 (2009) 
085501. 

T. J. Stasevich, Hailu Gebremariam, T.L. Einstein, M. 
Giesen, G. Steimer, and H. Ibach, Low-Temperature 
Orientation Dependence of Step Stiffness on {111} 
Surfaces, Phys. Rev. B 71 (2005) 245414. 

T. J. Stasevich, Ph. D. dissertation, Univ. of Maryland, 
2006 (unpublished) [http://drmn.lib.umd.edu/ 
bitstream/1903/4071/l/umi-umd-3818.pdf]. 

T. J. Stasevich, T. L. Einstein, R. K. P.Zia, M. Giesen, H. 
Ibach, and F. Szalma, The Effects of Next-Nearest- 
Neighbor Interactions on the Orientation Dependence 
of Step Stiffness: Reconciling Theory with Experiment 
for Cu(OOl), Phys. Rev. B 70 (2004) 245404. 

T. J. Stasevich, T. L. Einstein, and S. Stolbov, Ex¬ 
tended Lattice Gas Interactions of Cu on Cn(lll) and 
Gu(OOl): Ab-Initio Evaluation and Implications, Phys. 
Rev. B 73 (2006) 115426. 

J. Stewart, O. Pohland, and J. M. Gibson, Elastic- 
displacement field of an isolated surface step, Phys. 
Rev. B 49 (1994) 13848. 

1. N. Stranski, Zur Theorie des Kristallwachstums, Z. 
Phys. Chem. (Leipzig) 136 (1928) 259. 

W. P. Su, J. R. Schrieffer, and A. J. Heeger, Solitons 
in Polyacetylene, Phys. Rev. Lett. 42, 1698 (1979); 
Soliton excitations in polyacetylene, Phys. Rev. B 22 
(1980) 2099. 

S. Surnev, K. Arenhold, P. Coenen, B. Voigtlander, and 
H. P. Bonzel, Scanning tunneling microscopy of equi¬ 
librium crystal shapes, J. Vac. Sci. Technol. A 16 
(1998) 1059. 

B. Sutherland, C. N. Yang, and G. P. Yang, Exact Solu¬ 
tion of a Model of Two-Dimensional Ferroelectrics in 
an Arbitrary External Electric Field, Phys. Rev. Lett. 
19 (1967) 588. 

B. Sutherland, Quantum Many-Body Problem in One Di¬ 
mension: Ground State, J. Math. Phys. 12 (1971) 246; 
Exact Results for a Quantum Many-Body Problem in 
One Dimension, Phys. Rev. A 4 (1971) 2019. 

A. Szczepkowicz, A. Ciszewski, R. Bryl, G. Oleksy, C.- 
H. Nien, Q. Wu, and T. E. Madey, A comparison of 
adsorbate-induced faceting on flat and curved crystal 
surfaces. Surface Sci. 599 (2005) 55. 

K. Thiirmer, J.E. Reutt-Robey, and E. D. Williams, Nu- 
cleation limited crystal shape transformations. Surface 
Sci. 537 (2003) 123. 

E. Tosatti, private communication, March 2014. 

S. B. van Albada, M. J. Rost, and J.W. M. Frenken, 
Asymmetric and symmetric Wulff constructions of is¬ 
land shapes on a missing-row reconstructed surface, 
Phys. Rev. B 65 (2002) 205421. 

H. van Beijeren, Exactly Solvable Model for the Rough¬ 
ening Transition of a Crystal Surface Phys. Rev. Lett. 
38 (1977) 993. 



30 


H. van Beijeren and I Nolden, “The Roughening Transi¬ 
tion,” in: Structure and Dynamics of Surfaces II: Phe¬ 
nomena, Models, and Methods, edited by W. Schom- 
mers and P. von Blanckenhagen, Topics in Current 
Physics, vol. 43 (Springer, Berlin, 1987), pp. 259-300. 

J. Villain, “Two-Dimensional Solids and Their Interac¬ 
tions with Substrates” in: Ordering in Strongly Fluc¬ 
tuating Condensed Matter Systems, edited by T. Riste 
(Plenum, New York, 1980), pp. 221-260. 

M. von Laue, Der Wulffsche Satz fiir die Gleichgewichts- 
form von Kristallen, Z. Kristallogr., Mineral. Petrogr. 
105 (1943) 124. 

D. A. Walko and 1. K. Robinson, Energetics of oxygen- 
induced faceting on Cu(115), Phys. Rev. B 64 (2001) 
045412. 

Z. Wang and P. Wynblatt, The equilibrium form of pure 
gold crystals. Surface Sci. 398 (1998) 259. 

G. Wang, J. F. Webb, and J. Zi, The strictly attractive 
XjH? interaction between steps of crystal surfaces. Sur¬ 
face Sci. 601 (2007) 1944. 

X.-S. Wang, J. L. Goldberg, N. C. Bartelt, T. L. Einstein, 
and E. D. Williams, Terrace Width Distributions on 
Vicinal Si(lll), Phys. Rev. Lett. 65 (1990) 2430. 

J. D. Weeks, “The Roughening Transition” in: Ordering 
in Strongly Fluctuating Condensed Matter Systems, 
edited by T. Riste (Plenum, New York, 1980), pp. 293- 
317. 

J.D. Weeks, private discussions (2014). 

A. A. Wheeler, Cahn-Hoffman ^-Vector and Its Relation 
to Diffuse Interface Models of Phase Transitions, J. 
Stat. Phys. 95 (1999) 1245. 

E. D.Williams, R. J. Phaneuf, Jian Wei, N. C. Bartelt, 
and T. L. Einstein, Thermodynamics and statistical 
mechanics of the faceting of stepped Si(lll), Surface 
Sci. 294 (1993) 219; erratum 310 (1994) 451. 

E. D. Williams and N. C. Bartelt, “Thermodynamics and 
Statistical Mechanics of Surfaces” in: Physical Struc¬ 
ture of Solid Surfaces, edited by W. N. Unertl (Else¬ 
vier, Amsterdam, 1996), Handbook of Surface Science, 
vol. 1, S. Holloway and N.V. Richardson, series eds., 
pp. 51-99. 

E. D. Williams and N. C. Bartelt, Surface Faceting and 


the Equilibrium Crystal Shape, Ultramicroscopy 31 
(1989) 36. 

P. E. Wolf, F. Gallet, S. Balibar, E. Rolley, and P. 
Nozieres, Crystal growth and crystal curvature near 
roughening transitions in hep ^He, J. Phys. (France) 
46 (1985) 1987. 

P. E. Wolf, S. Balibar, and F. Gallet, Experimental Ob¬ 
servation of a Third Roughening Transition on hep "^He 
Crystals, Phys. Rev. Lett. 51 (1983) 1366. 

M. Wortis, “Equilibrium Crystal Shapes and Interfacial 
Phase Transitions,” in: Chemistry and Physics of Solid 
Surfaces, VH, edited by R. Vanselow and R. Howe 
(Springer-Verlag, Berlin, 1988), pp. 367-405. 

G. Wulff, Zur Erage der Geschwindigkeit des Wachstums 
und der Auflosung der Krystallflachen, Z. Krystallo- 
graphie und Mineralogie 34 (1901) 449. 

M. Yamamoto, K. Sudoh, H. Iwasaki, and E. D. Williams, 
Anomalous decay of multilayer holes on SrTiO3(001), 
Phys. Rev. B 82 (2010) 115436. 

C. P. Yang, Exact Solution of a Model of Two- 
Dimensional Eerroelectrics in an Arbitrary External 
Electric Eield, Phys. Rev. Lett. 19 (1967) 586. 

M. Yoon, S. G. J. Mochrie, D. M. Zehner, G. M. Watson, 
and D. Gibbs, Eaceting and the orientational phase 
diagram of stepped Pt(OOl) surfaces, Phys. Rev. B 49 
(1994) 16702. 

K. Yosida, Magnetic Properties of Gu-Mn Alloys, Phys. 
Rev. 106 (1957) 893. 

D. K. Yu, H. P. Bonzel, M. SchefFler, The stability of vic¬ 
inal surfaces and the equilibrium crystal shape of Pb 
by first principles theory. New J. Phys. 8 (2006) 65. 

R. K. P. Zia, “Anisotropic surface tension and equilibrium 
crystal shapes,” in: Progress in Statistical Mechanics, 
edited by G. K. Hu (World Scientific, Singapore, 1988), 
pp. 303-357. 

R. K. P. Zia, Exact equilibrium shapes of Ising crystals 
on triangular/honeycomb lattices, J. Stat. Phys. 45 
(1986) 801. 

J.-K. Zuo, T. Zhang, J. E. Wendelken, and D. M. Zehner, 
Step bunching on TaC(910) due to attractive step-step 
interactions, Phys. Rev. B 63 (2001) 033404. 


[1] This term was coined by a speaker at a workshop in Tra¬ 
verse City in Angust 1996—see Duxbury & Pence [1997] 
for the proceedings—and then used by several other speak¬ 
ers. 

[2] As exposited clearly in Callen [1985], one considers a 
[convex] function y = y{x) and denotes its derivative as 
p = dyjdx. If one then tries to consider p instead of 
X as the independent variable, there is information lost: 
one cannot reconstruct y(x) uniquely from y(p). Indeed, 
y = y{p) is a- first-order differential equation, whose in¬ 
tegration gives y = y(x) only to within an undetermined 
integration constant. Thus, y ~ y{p) corresponds to a fam¬ 
ily of displaced curves, only one of which is the original 
y = y{x). The key concept is that the locus of points satis¬ 


fying y = y(x) can be equally well represented by a family 
of lines tangent to y{x) at all x, each with a ^-intercept tp 
determined by the slope p at {x, y{x)). That is, tp = tp{p) 
contains all the information of y = y{x). Recognizing that 
p = {y — Ip)/{x — 0), one finds the transform tp = y — px. 
Readers should recall that this is the form of the relation¬ 
ship between thermodynamic functions, particularly the 
Helmholtz and the Gibbs free energies. 

[3] The conventional designation of this exponent is A or 0. 
However, these Greek letters are the Lagrange multiplier 
of the ECS and the polar angle, respectively. Hence, we 
choose d for this exponent. 

[4] There are some minor differences in prefactors from 
Bartelt et al. [1989] 



